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ABSTRACT 


This  report  covers  work  of  the  Optics  Division  at  Lincoln  Laboratory 
for  the  period  1  January  through  30  June  1972.  The  topics  covered  are 
laser  technology  and  propagation,  optical  measurements  and  instrumen-* 
tation,  and  laser  radar. 

Additional  information  on  the  optics  program  may  be  found  in  the  semi¬ 
annual  technical  summary  reports  to  the  Advanced  Research  Projects 
Agency. 


Accepted  for  the  Air  Force 
Joseph  J.  Whelan,  USAF 

Acting  Chief,  Lincoln  Laboratory  Liaison  Office 


CONTENTS 


Abstract  iii 

Introduction  vii 

Reports  on  Optics  Research  xi 

Organization  xv 

L  LASER  TECHNOLOGY  AND  PROPAGATION  1 

A.  CW  Theory  1 

1.  Laser  Propagation  Through  a  Turbulent  Medium  1 

2.  Interaction  of  Heating  and  Turbulence  in  Thermal  Blooming  10 

3.  Calculations  of  Phase  Correction  for  a  Laboratory  Thermal 

Blooming  Experiment  13 

B.  Effects  15 

1.  Laser-Induced  Air  Breakdown  15 

2.  Surface  Interactions  19 

3.  Explosive  Vaporization  of  Fog  Droplets  by  a  10.6-fxm 

Laser  Pulse  31 

C.  Devices  37 

1.  Electron  Beam  (e-Beam)  37 

2.  The  Possibility  of  Generating  Laser  Plasmas 

by  Photoionization  51 

3.  Other  Device  Work  54 

II.  OPTICAL  MEASUREMENTS  AND  INSTRUMENTATION  59 

A.  ZnS  Up-Conversion  59 

B.  IR  Camera  (4  inch)  '  59 

C.  Interferometric  Imaging  Techniques  62 

1.  Introduction  62 

2.  Experimental  Program  63 

3.  Results  65 

4.  Conclusions  67 

D.  Airborne  Image  Studies  68 

1,  Measurements  of  the  Tropopause  Effect  68 

2,  Boundary  Layer  71 

E.  Tunable  Infrared  Laser  Spectroscopy  of  Sulfur  Dioxide  72 

1 .  Introduction  72 

2.  Diode  Laser  Fabrication  and  Experimental  Procedure  73 

3.  Theory  73 

4.  Results  74 

5.  Application  of  Theoretical  Spectra  to  Pollutant  Gas 

Monitoring  77 

6.  Acknowledgments  78 

F.  Extension  of  Tunable  Laser  Techniques  to  Air  Pollution 

Detection  78 

1.  ,Spectrophone  Investigation  79 

2.  In -Situ  Auto  Exhaust  Measurements  79 

3.  Long-Path  Laser  Transmission  Studies  83 


V 


Contents 


III.  LASER  RADAR  AND  TRACKING 

A.  Instrument  for  Atmospheric  Modulation  Transfer  Function 
Measurements 

1 ,  Introduction 

2.  MTF  Electronics 

B.  Automatic  Range  Tracker 

C.  Atmospheric  Optical  Coherence  Measurements 
at  10.6  and  0.63  pm 

D.  InSb  Isolator 


85 

85 

85 

86 

87 

88 
89 


y 


VI 


INTRODUCTION 


I.  LASER  TECHNOLOGY  AND  PROPAGATION 

The  computer  code  for  calculating  thermal  blooming  described  previously  has  been  extended 
to  include  the  effects  of  turbulence  in  the  atmosphere. 

A  theoretical  study  has  been  made  of  two  mechanisms  which  have  been  suggested  as  impor¬ 
tant  in  thermal  blooming  calculations:  the  dependence  of  a  upon  localized  temperature  in  a 
thermally  turbulent  medium,  and  the  effect  of  velocity  turbulence  which  would  tend  to  cool 
off  the  heated  air.  The  study  shows  that  for  power  levels  presently  contemplated,  both  these 
effects  are  negligible. 

A  series  of  propagation  calculations  have  been  performed  which  show  that  CW  thermal  bloom¬ 
ing  can  be  considerably  reduced  by  an  initial  phase  correction  of  the  laser  beam. 

Several  experiments  were  performed  to  investigate  the  effects  of  atmospheric  dust  on  the 
breakdown  threshold  of  10.  b-pm  laser  beams  with  large  focal  spots.  The  lasers  used  in  these 
experiments  were  a  1-kJ  pulsed  e-beam  at  Avco  Everett  Research  Laboratories  (AERL),  and 
a  3-J  pulsed  TEA  laser  at  Lincoln  Laboratory.  The  main  conclusions  from  these  experiments 
are:  the  breakdown  threshold  of  clean  air  is  close  to  that  predicted  by  microwave  breakdown 
theory;  dust  particles  can  lower  the  threshold  below  this  value  roughly  in  proportion  to  their 
size;  and  the  presence  of  atmospheric  dust  can  explain  the  ’’spot  size  effect”  reported  recently 
by  D.  C.  Smith, 

Impulse  data  obtained  with  the  AERL  1-kJ  pulsed  CO^  laser  have  been  roughly  correlated  with 
a  simple  cylindrical  blast  wave  model.  The  apparent  saturation  of  impulse  with  increasing 
laser  energy  was  observed  in  the  presence  of  multiple  aerosol-induced  breakdowns  in  the 
laser  beam. 

Observations  of  plasmas  produced  by  irradiation  of  surfaces  with  the  AERL  pulsed  CO2  laser 

were  carried  out.  The  plasma  and  associated  shock-wave  expansions  were  observed  with  a 

6  ° 
framing  camera  (2  x  10  frames/sec).  The  plasma  absorption  at  wavelengths  4881 A  and 

10.  6  pm  was  measured  and  used  to  infer  electron  density  histories. 

The  explosive  vaporization  of  individual  water  droplets  (5- to  25-|JLm  radius)  by  a  10.  b-fj-m 
laser  pulse  has  been  observed  with  a  high-speed  schlieren  photography  system.  The  hot  vapor 
sphere  and  the  shock  wave  produced  by  the  explosion  can  be  clearly  distinguished  in  the  schlie¬ 
ren  photographs. 

The  Lincoln  500-]  CO 2  e-beam  excited  laser  is  in  the  final  stages  of  construction.  Laser  op¬ 
erations  should  commence  in  August  1972,  and  full  device  operation  is  expected  during  October. 

The  rate  equations  and  temperature  balance  equation  have  been  solved  to  predict  small  signal 
gain,  pulse  energy,  and  pulse  length  for  an  e-beam  excited  CO2  laser.  The  initial  spike,  as 
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observed  experimentally  for  the  laser  pulse  from  the  Avco  e-beam  device,  has  been  predicted 
but  the  computed  temporal  width  is  narrower  than  that  observed. 

A  2.  5-m  unstable  confocal  positive  branch  resonator  has  been  designed  for  the  pulsed  500-J 
CO^  laser  so  that  effective  use  can  be  made  of  the  large-volume  short-length  (Im)  gain  me¬ 
dium.  Although  several  different  outcouplings  and  geometries  will  be  tried,  it  is  expected  that 
the  optimum  mode  quality  outcoupling  85  percent  will  produce  45  percent  of  the  net  energy 
within  the  first  Airy  dark  ring. 

The  Freedholm  integral  equation  is  used  to  determine  the  asymmetric  unstable  resonator  cav-  t 

ity  mode  field  distribution  without  a  gain  medium.  Conversion  is  made  by  function  expansion 
to  the  eigenvalue  solution  of  an  infinite  matrix.  A  finite  matrix  solution  approximation  is  then 
used  for  computer  calculation. 

A  3-nsec  600-keV  e-beam  source  has  been  used  to  initiate  a  CO2  laser  beam.  Output  optical 
pulses  up  to  120fisec  and  22.  4J  have  been  observed  from  an  active  volume  of  about  1  liter. 

A  pulsed  CO2  laser  beam  has  been  shaped  in  time  by  passing  the  beam  through  an  absorption 
cell  that  is  external  to  the  laser  and  filled  with  SF^  or  Freon-12.  The  initial  spike  on  the 
laser  pulse  has  been  successfully  attenuated  when  passed  through  the  absorption  cell. 

A  theoretical  investigation  has  been  carried  out  to  determine  the  difficulties  of  amplifying  very 
short  ('-I  nsec)  pulses  or  frequency-chirped  pulses  at  10-|jLm  optical  wavelength  to  high  ener¬ 
gies  (~200J/sequence).  The  conclusions  indicated  that  only  the  volume  preionized  devices  are 
reasonable  in  the  last  stages  of  amplification.  Considerations  are  made  of  the  saturation  and 
pulse  sharpening  limitations  on  such  devices.  Frequency  distortion  and  index  changes  are  not 
seen  to  be  a  serious  problem.  The  most  severe  difficulties  are  found  to  be  in  bulk  medium 
heating  and  cathode  shock-wave  distortion  as  limiting  the  coherence  time  of  the  pulses. 

12  3 

It  should  be  possible  to  produce  electron  densities  on  the  order  of  10  /cm  at  atmospheric 
pressure  by  photoionization  in  a  laser  gas  such  as  CO2  or  CO.  The  ionization  produced  will  be 
uniform  throughout  large  laser  medium  volumes  if  the  photon  absorption  length  is  large,  and 
this  can  be  accomplished  by  seeding  the  laser  gas  with  a  low  partial  pressure  of  gas  that  has  a 
low  ionization  potential.  This  will  then  allow  deep  photon  penetration  and  low-energy  UV  pho¬ 
tons  to  be  utilized  for  ionization. 

A  small  ionization  cell  containing  a  Xe  flash  lamp  with  suprasil  window  and  two  flat  parallel 

electrodes  is  used  to  study  volume  photoelectron  ionization.  When  0.2torr  of  tripropylamine 

10  3 

is  added  to  500-torr  N2,  an  electron  density  n^  =  10  /cm  is  measured  upon  flash-lamp  igni¬ 
tion.  The  current  is  found  to  vary  linearly  with  applied  electrode  voltage,  which  implies  an  n^ 
production  only  by  photoionization. 

The  small  signal  gain  and  laser  level  pumping  rates  for  a  CO2  laser  have  been  scaled  with 
pressure  to  10  atm.  Using  longitudinal  excitation  by  a  10^-A,  3-nsec  e-beam  source  (Febatron 

706),  the  peak  gain  for  the  P2Q  line  has  been  calculated  to  be  2.  2  percent/cm.  » 

A  20-cm^  CO2  laser  device  has  been  constructed  to  operate  at  10 atm,  and  experiments  with 
this  laser  will  commence  when  the  Febatron  706  is  available  as  a  longitudinal  e-beam  source. 
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Introduction 


Measurements  and  evaluation  of  data  have  been  completed  on  the  Lamb  dip  for  16  different 
rotational-vibrational  transitions  in  sealed-off  CO  lasers  operating  at  a  total  pressure  of 

l.Storr. 


II.  OPTICAL  MEASUREMENTS  AND  INSTRUMENTATION 

In  the  ZnS  upconversion  process,  AC  effects  have  been  measured  using  chopped  IR  radiation 
as  a  source  and  constant  UV  illumination  for  excitation.  The  AC  SNR  is  1  to  2  orders  of  mag¬ 
nitude  worse  than  the  SNR  obtained  in  previous  measurements  where  the  phosphor  is  exposed 
to  DC  IR  radiation  after  (not  during)  UV  excitation.  Apparatus  to  determine  the  S/N  charac¬ 
teristics  of  the  phosphor  screen  resolution  elements  is  being  constructed. 

Measurements  have  been  made  relevant  to  the  imaging  quality  of  a  cooled  4 -inch -diameter, 
3-element  refractor  camera.  Tests  have  demonstrated  the  mechanical  integrity  of  the  design 
and  its  diffraction-limited  performance  in  the  wavelength  range  of  7.5  to  12(JLm.  Extraneous 
radiation  at  angles  about  18°  off  the  optical  axis  has  been  observed,  and  some  effort  has  been 
made  to  determine  the  cause  of  these  peaks. 

The  shearing  interferometer  was  employed  on  board  the  NKC-135  aircraft  with  auxiliary  op¬ 
tics  suspended  from  a  window-mounted  air-foil  plate  to  measure  the  modulation  transfer  func¬ 
tion  of  the  local  boundary  layer  within  25  cm  of  the  aircraft  skin.  The  existence  of  image  deg¬ 
radation  due  to  a  boundary  layer  was  established. 

Stellar  images  recorded  in  high-altitude  flight  with  a  diffraction-limited,  long-focal-length 
telescope  have  been  analyzed  to  determine  the  possible  effects  of  the  tropopause  in  image  deg¬ 
radation.  The  results  suggest  that  the  difference  between  the  stratosphere  and  troposphere 
is  not  measurable  in  the  presence  of  boundary -layer  turbulence,  and  a  determination  of  turbu¬ 
lence  effects  agrees  quite  well  with  interferometer  results  measured  separately. 

Detailed  tunable  diode  laser  studies  have  been  made  of  the  band  of  SO^  centered  around 
8.7|JLm.  Comparison  between  the  experimental  results  and  theoretical  predictions  based  on 
microwave  data  has  resulted  in  the  establishment  of  a  more  accurate  value  for  the  band  center 
frequency.  Characteristics  of  the  absorption  spectra  under  varying  conditions  of  temperature 
and  pressure  are  used  to  determine  optimum  wavelengths  for  monitoring  applications. 

Several  extensions  of  tunable  diode  laser  techniques  to  air-pollution  monitoring  are  described. 
The  sensitivity  of  acousto-optic  absorption  cells  and  direct  transmission  measurements  for 
air-pollution  monitoring  have  been  compared;  real-time,  quantitative,  in -situ  measurements 
of  auto  exhaust  have  been  made  and  are  reported;  and  long  path  laser  transmission  studies  of 
ambient  air  pollution  are  described. 


III.  LASER  IMAGING  RADAR 

The  10.  6-^im  imaging  radar  has  been  used  to  measure  cross-section  characteristics  of  two 
satellite  models  and  a  Terrier  radome. 
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Introduction 


A  simple,  portable  modulation-transfer-function  (MTF)  instrument  based  on  a  Sieman's  star 
reticle  has  been  built  and  tested.  Under  laboratory  conditions,  the  instrument  has  demon¬ 
strated  diffraction-limited  performance.  In  the  field,  the  atmospheric  MTF  has  been  meas¬ 
ured  under  a  wide  variety  of  conditions  in  support  of  atmospheric  spatial  coherence  measure¬ 
ments  which  were  being  made  with  a  48-inch  telescope  at  10. 6 pm. 

Electronics  have  been  constructed  both  for  control  of  the  MTF  instrument  and  for  processing 
of  the  output  signal.  The  MTF  processor  provides  a  calibrated  output  voltage  proportional  to 

the  percentage  modulation.  Various  control  circuits  for  the  Firepond  laser  radar  have  also  » 

been  designed  and  implemented. 

Spatial  coherence  measurements  have  been  made  at  10.6  and  0.63|jLm.  The  data  at  10.6  pm 
are  currently  being  reduced,  but  qualitatively  agree  with  the  results  obtained  at  0.63pm. 

These  latter  data,  taken  over  a  5.4-km  horizontal  path  at  an  average  68m  above  the  ground, 
indicate  that  the  coherence  diameter  r  (0. 63  pm)  for  a  wide  variety  of  atmospheric  conditions 
varies  from  1cm  <r^  (0.63pm)  <20cm.  From  this,  we  predict  that  the  10.  6-pm  data  will 
reveal  30cm  <  r^(10.6pm)  <  6m  for  this  path. 

A  10. 6-pm  isolator  has  been  fabricated  using  Faraday  free-carrier  rotation.  The  isolator  has 
a  forward  insertion  loss  of  2dB,  a  total  reverse  loss  of  22  dB,  and  a  maximum  throughput 
power  of  30  W  CW. 
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L  LASER  TECHNOLOGY  AND  PROPAGATION 


A.  CW  THEORY 

1.  Laser  Propagation  Through  a  Turbulent  Medium 

In  previous  reports,  we  discussed  the  propagation  of  a  high-power  laser  beam  through  a 
nonturbulent  atmosphere.  In  the  past  six  months,  we  have  extended  our  capabilities  to  include 
propagation  through  a  turbulent  atmosphere. 

The  standard  approach  to  the  calculation  of  propagation  of  a  low- power  beam  through  a  tur¬ 
bulent  medium  is  by  analytic  computation  of  ensemble  average  values  of  various  quantities.  (A 
good  review  of  this  field  is  found  in  Ref.  1.)  However,  since  analytic  methods  had  not  been  suc¬ 
cessful  in  the  high-power,  nonturbulent  case,  another  approach  was  needed.  We  accomplished 
this  by  simulating  individual  realizations  of  the  atmosphere,  with  the  constraint  that  the  ensem¬ 
ble  of  these  realizations  should  have  the  statistical  properties  of  an  accepted  (or  at  least,  widely 
used)  model  of  the  turbulent  atmosphere.  In  the  following  section,  we  argue  that  the  velocity 
turbulence  of  the  atmosphere  is  of  little  importance  except  in  the  case  of  an  extremely  high- 
power  beam,  and  that  for  the  majority  of  cases  it  is  an  adequate  approximation  simply  to  add  the 
naturally  occurring  temperature  (or  index)  fluctuations  to  the  temperature  (or  index)  change  in¬ 
duced  by  the  beam  itself.  Since  our  numerical  method  is  capable  of  computing  the  propagation 
through  an  arbitrary  atmosphere  (so  long  as  the  index  gradients  are  not  excessive),  the  only 
remaining  problem  is  the  generation  of  suitable  realizations  of  the  temperature  fluctuations. 

For  our  purposes,  sufficient  information  about  the  statistical  properties  of  the  index  of  re¬ 
fraction  is  contained  in  the  autocorrelation  function  <(n(  r )  n(  r  +  p  )]>.  If  the  atmosphere  is  con¬ 
sidered  to  be  statistically  homogeneous,  this  function  is  independent  of  r,  and  can  be  described 
equally  well  by  the  (three-dimensional)  spectrum 
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(Ztt) 
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This  spectrum  is  generally  used  for  the  statistical  description  of  the  atmosphere.  But  in  our 
propagation  calculations,  the  quantity  n  does  not  enter  directly,  but  only  as  an  integral,  J ndz. 
To  this  integral  corresponds  a  two-dimensional  spectrum: 
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It  is  easy  to  show  that  ^2  related  to  by  the  equation 
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where  jQ(u)  e  sinu/u. 
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Fig.  1-1 .  Iso-irrad lance  contour  plots  for  conditions  corresponding  to  FIrepond  facility:  122-cm 
Cassegrainlan  (0,33  obscuration  ratio)  focused  at  5.4  km,  10.6  pm.  Turbulence  Is  assumed 

2  -14  -2/3 

to  be  homogeneous  and  isotropic,  with  C  =2X10  m  ^  .  Four  different  real  Izatlons.  Grid 
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We  now  wish  to  generate  values  of  J ndz  on  a  finite  two-dimensional  grid.  We  do  this  by 
generating  an  array  of  random  complex  numbers  ^  with  zero  mean  and  Gaussian  distribution 
in  the  space  k  ,  k  multiplying  by  the  required  function  of  k  k  and  taking  a  Fourier  transform: 

X  y  ^  y 
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where  and  are  the  total  widths  of  the  grid  in  the  x-  and  y-directions,  respectively. 

Requiring  a  least- mean- square  fit  of  J  ndz  based  on  a  finite  grid  to  the  continuum  expres¬ 
sion,  one  can  show  that  a  convolution  is  necessary  in  order  to  take  proper  account  of  frequencies 
present  in  the  continuous  spectrum  t)ut  not  present  in  the  spectrum  for  the  discrete  wave- 


numbers  k  =  27rm/L  and  k  =  27rn/L  . 
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We  use  the  scheme  described  above  to  generate  an  array  of  phase  changes  (proportional  to 
J  ndz)  twice  for  every  step  in  z,  in  accordance  with  our  method  of  calculation.^  The  arrays  are 
taken  to  be  statistically  independent  of  one  another,  since  the  z-steps  actually  used  are  consider¬ 
ably  larger  than  the  outer  scale  of  turbulence,  except  in  the  immediate  vicinity  of  a  focus  where 
the  small  optical  lever  arm  makes  phase  changes  of  little  importance.  For  the  same  reason, 
we  believe  that  the  assumption  of  a  Gaussian  distribution  for  the  phase  changes  is  justified.  For 
the  spectrum  we  have  used  the  modified  von  Karman  spectrum  in  the  form 


^^(k)  =  0.033 
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where  L  is  the  outer  scale,  and  is  the  inner  scale  of  turbulence, 
o  o 

Figure  I-l  shows  some  iso-irradiance  contour  plots  that  result  when  we  apply  our  method 
to  a  particular  experimental  situation,  without  nonlinear  effects.  Each  plot  corresponds  to  an 
instantaneous  picture  of  the  beam,  and  the  variation  from  one  plot  to  another  is  due  to  the 
statistical  fluctuations  of  the  assumed  atmosphere. 

As  overall  measures  of  the  effect  of  turbulence,  we  have  chosen  two  quantities.  One  is  the 
equivalent  area,  defined  as 

A  =  ^/ 

2  /  dxdy 

For  well-behaved  intensity  patterns,  the  equivalent  area  so  defined  agrees  with  the  area  defined  by 

■^rms  ^  J  [{x  -  X  -Ky  -  y  )^]  Idxdy/^  Idxdy 

based  on  the  standard  deviations.  But  for  the  Airy  pattern,  the  rms  area  diverges  whereas  the 
equivalent  area  gives  a  reasonable  value.  In  addition,  the  standard  deviations  cannot  be  reliably 

3 

calculated  for  noisy  experimental  data,  since  the  weighting  of  the  intensity  with  r  (as  can  be 
seen  by  putting  the  integral  into  cylindrical  coordinates)  leads  easily  to  numerical  difficulties 
for  large  r.  We  therefore  prefer  to  use  the  equivalent  area.  This  is  usually  divided  by  the 
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corresponding  quantity  in  the  absence  of  turbulence  and  a  nondimensional  "beam  spread"  is 
defined: 

cr  ?  =  (A  “A  )/A 
A  '  eq  o  '  o 

The  other  quantity  is  the  displacement  of  the  center  of  irradiance  of  the  pattern: 
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and  the  nondimensional  "wander"  (squared)  is  defined  as 

(7  ^  =  7r(w  ^  +  w  ^  )/A 
w  ~  X  y  o 

Z 

In  Figs.  1-2  and  1-3,  we  plot  the  values  of  these  quantities  as  a  function  of  C  the  measure  of 

2  ^ 

the  strength  of  turbulence.  At  a  given  value  of  C^,  the  different  points  correspond  to  different 
realizations  of  the  atmosphere.  Points  that  are  connected  by  dashed  lines  correspond  to  reali- 
zations  that  differ  only  by  the  value  of  C^;  i.e.,  the  random  numbers  used  to  generate  the  reali¬ 
zations  are  the  same. 

The  fact  that  these  dashed  lines  are  nearly  parallel  to  one  another  suggests  that  it  may  be 
possible  to  study  the  dependence  of  turbulent  propagation  on  various  parameters  by  maintaining 
a  single  set  of  random  numbers  and  varying  the  parameters,  without  having  to  accumulate  a 
large  number  of  statistics  for  each  set  of  parameters.  This  possibility  is  made  more  plausible 
by  an  examination  of  the  contour  plots  of  Figs.  1-4  and  1-5,  which  result  from  just  such  a  study, 
where  a  set  of  dimensionless  numbers  appropriate  for  this  problem  have  been  chosen  and  varied 
systematically.  The  dimensionless  numbers  are  the  Fresnel  number 


N  E  k  a  /L 
F  o  m' 


the  f /number 

f/  E  L/2a 
^  '  m 


the  "turbulent  distortion  number" 

N„=a  (k^C^L)^/® 

T  m  o  n 

and  the  ratios  a  /L  and  a  /f  .  Here,  a  is  a  characteristic  radius  in  the  initial  plane,  k  is 
o  m^  o  m  "  o 

the  optical  wavenumber,  and  both  the  initial  radius  of  curvature  and  the  distance  to  the  plane  of 
observation  are  taken  to  be  L.  For  Figs.  1-4  and  1-5,  the  parameter  has  been  kept  constant, 
while  the  other  parameters  are  varied.  From  the  "family  resemblance"  of  the  plots,  it  seems 
clear  that  scaling  laws  can  be  deduced  in  this  way  and  that  is  the  most  important  parameter 
in  the  range  studied.  Its  effect  on  the  beam  spread  and  wander  are  already  visible  in  Figs.  1-2 
and  1-3;  Figs.  1-6  and  1-7  show  the  effects  of  varying  Np. 

We  can  now  proceed  to  compute  the  propagation  of  a  high-power  beam  through  a  turbulent 
atmosphere.  In  Fig.  1-8,  we  see  a  comparison  of  results  for  a  high-power  beam  propagating 
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Fig.  1-8.  Thermal  blooming  in  presence  of  turbulence.  Top  contour  plot  shows  beam  with  turbulence 
absent  and  thermal  blooming  present.  Left  three  contour  plots  show  beam  in  absence  of  thermal  bloom¬ 
ing,  after  propagation  through  three  realizations  of  turbulent  atmosphere  (C^  =  10“^^  m  “2/3).  Right 
three  contour  plots  show  thermally  blooming  beam  after  passing  through  same  three  realizations  of  tur¬ 
bulent  atmosphere.  Initial  beam  was  an  infinite  Gaussian,  l/e^  power  diameter  1  m,  focused  and  ob¬ 
served  at  range  of  10  km,  X=  10.6  pm,  and  aP/v  =  10“4. 

w 
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Fig.  I“9.  Total  loss  (thermal  blooming  plus  turbulence)  as  a  function  of  slewing  rate-  Solid  line  shows 

loss  In  absence  of  turbulence,  while  hatched  areas  show  range  of  loss  associated  with  various  values  of 

C^,  Infinite  Gaussian,  l/e^  power  diameter  1  m,  range  3  km,  X  =  20  pm,  aP/v  =  10"'^. 
n  w 

through  the  same  turbulent  atmosphere.  Figure  1-9  shows  the  intensity  loss  when  turbulence 
and  thermal  blooming  are  combined,  for  a  particular  example.  More  detailed  studies  are 

under  way.  C.  Bradley 

J.  Herrmann 

2.  Interaction  of  Heating  and  Turbulence  in  Thermal  Blooming 

In  our  calculation  of  the  propagation  of  a  high-power  beam  through  a  turbulent  atmosphere, 
we  assumed  that  the  temperature  variations  due  to  heating  may  be  simply  added  to  those  already 
present  in  the  atmosphere.  It  has  been  suggested  that  the  interaction  between  heating  and  tur¬ 
bulence  may  give  rise  to  additional  temperature  variations  sufficiently  large  that  the  propagation 
will  be  significantly  affected.  Two  mechanisms  have  been  proposed: 

(a)  Warmer  blobs  have  a  larger  absorption  coefficient  than  cooler  ones, 
so  heating  will  tend  to  amplify  the  variations. 

(b)  Velocity  fluctuations  in  the  medium  will  tend  to  cool  off  the  heated 

air  more  effectively,  and  also  will  cause  further  temperature  fluctuations. 

Our  purpose  here  is  to  estimate  these  effects  and  show  that  both  are  negligible,  at  least  at  pres¬ 
ently  contemplated  power  levels. 

We  assume  that  wind  velocities  are  much  smaller  than  the  sound  velocity,  so  that  the  pres¬ 
sure  may  be  considered  constant.  The  temperature  is  thus  considered  a  passive  additive  to  the 
turbulent  velocity  field.  If  the  winds  are  not  so  slow,  the  heating  effects  will  be  negligible  any¬ 
how.  Then,  the  equation  of  energy  conservation  may  be  written  as 
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VT  = 
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(I-l) 


We  assume  the  properties  of  the  radiation  source  to  be  independent  of  time  (or  slowly  changing), 
and  divide  the  properties  of  the  medium  into  mean  and  fluctuating  parts: 


T  =  T  +  T . 
o  1 


V  =  V  +  V . 

o  1 
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The  mean  parts  are  assumed  to  satisfy 
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To  this  might  be  added  a  term  ~v.  •  VT^,  but  we  argue  below  that  this  term  is  negligibly  small 
for  practical  application.  Then  we  have 
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(1-3) 


We  now  separate  that  part  of  the  random  temperature  field  which  can  be  considered  as  the  frozen 
turbulence  T^,  satisfying  the  Taylor  h3rpothesis 


8T 


Then,  with  T^  =  we  get 
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In  our  propagation  code  (see  Sec.  A-1 )  in  a  turbulent  atmosphere  the  contributions  of  T^  and  T^ 
have  been  included.  We  must  now  estimate  the  contributions  on  the  right-hand  side  of  Eq.  (1-4) 
to  show  that  T^  is  negligible. 

The  time  for  a  fluid  particle  to  cross  the  beam  is  ~a/v^,  where  a  is  the  beam  radius,  so 
the  increment  in  T^  is 


AT., 
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We  may  compare  the  increment  from  the  first  term  (A^T^)  to  the  rms  value  of  T^  in  the  unper¬ 
turbed  fluid,  which  we  call  T  : 

f  r 


Now, 


A  T 

1^2  J_  ^  JL 

T  ~  T  V  D  c 
r  r  o  p 


(-) 

Lo 

(-) 

P  o 


E  _ r 

kT  T 


11 


Section  I 


where  E  is  the  lower  level  of  the  laser  transition. 


and  I  P/7ra 


so  we  find 


^^2 
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For  fairly  extreme  cases,  this  may  amount  to  percent,  which  we  consider  to  be  negligible. 

It  is  not  appropriate  to  compare  the  increment  from  the  second  term  {A^T^)  with  T^,  since 
under  conditions  of  low  lapse  rate  (at  dawn  or  dusk)  it  is  possible  to  have  large  velocity  fluctua¬ 
tions  and  small  temperature  fluctuations.  Instead,  we  compare  it  with  AT^,  the  change  in 
in  traversing  the  beam.  It  is  easily  seen  that  (A^T^/AT^)  ~  (v^/v^),  where  is  the  rms  value 

of  V,.  This  is  an  overestimate  of  the  effect,  however,  for  the  important  quantity  for  propagation 
1  .  . 
is  not  the  change  in  temperature  itself,  but  its  integral  over  the  length  of  the  beam,  which  is 

proportional  to  the  phase  change.  Since  v^  is  uncorrelated  over  paths  greater  than  some  outer 

scale  L  ,  while  T  adds  coherently  over  the  whole  path,  which  we  assume  is  much  larger  than 
0  0 

L^,  we  have 

r  AT  dz  cc  (AT  )  (Az) 

Jo  o 

J  A^T^dz  c  (A2T3)  (L^Az)1/2 


and  the  ratio  of  the  corresponding  phase  changes  is 

hi 

^^o  ~  ""o 

At  most,  (v  /v^)  ~  0.1,  while  (L^/Az)  ~  0.01  for  the  cases  in  which  we  are  interested.  Thus, 
we  consider  that  this  term  too  can  be  neglected. 

The  third  term  (A^T^)  describes  the  breakup  of  frozen  turbulence.  Since  VT^  ~  (T^/L^), 
we  find 


which  is  negligibly  small. 

The  fourth  term  (A^T^)  is  the  most  difficult  with  which  to  deal.  The  time  variation  of  T^  is 
due  to  the  fact  that  turbulence  in  the  initial  parts  of  the  beam  leads  to  scintillations  (time 
changes  in  I)  at  later  points  in  the  beam.  In  general,  the  rapidly  varying  part  of  I  has  high 
spatial  frequencies  and  therefore  tends  to  integrate  out  in  the  equation  for  T^,  whereas  the  low 
spatial  frequencies,  which  are  the  most  important  in  contributing  to  T^,  are  slowly  varying. 
Furthermore,  in  a  focused  beam,  the  relevant  time-scale  a/v^  becomes  smaller  as  the  beam 
travels  farther  and  as  the  scintillations  thus  become  larger.  Therefore,  A^T^  may  not  be  as 
important  as  it  appears  at  first  glance.  In  an  attempt  to  make  these  statements  more  quantitative, 
we  are  now  carrying  out  ’’computer  experiments”  on  low-power  beams. 

The  above  considerations  apply  to  a  beam  of  moderate  dimensions  (<1  m)  propagating  through 
the  atmosphere.  Recently,  Gebhardt,  et  have  studied  thermal  blooming  in  a  turbulent  me¬ 
dium  in  the  laboratory.  They  have  attempted  to  .make  v^  and  T^  as  small  as  possible,  and  the 
integral  scale  for  their  turbulent  medium  is  on  the  order  of  a  beam  diameter.  For  such  an 
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environment,  it  is  no  longer  permissible  to  neglect  the  term  -v^  •  VT^  in  Eq.  (1-2);  the  -us-ual 
practice  is  to  replace  it  by  an  eddy  diffusion  term  V  •  ^  While  this  is  no 

doubt  appropriate  for  their  experiment,  we  believe  that  our  order-of-magnitude  analysis  above 
shows  that  this  term  will  not  be  important  in  the  open  atmosphere. 

J.  Herrmann 
L.  C.  Bradley 

3.  Calculations  of  Phase  Correction  for  a  Laboratory  Thermal 
Blooming  Experiment 

In  order  to  design  a  laboratory  experiment  for  the  study  of  the  elimination  of  thermal  bloom¬ 
ing  by  an  initial  phase  correction,  a  series  of  propagation  calculations  has  been  performed. 

Figure  I-IO  shows  the  radius  of  the  unbloomed  beam  near  the  focus  for  six  aperture  sizes. 
All  cases  assume  a  Gaussian  beam  truncated  at  this  aperture  size  with  a  radius  of  curvature  of 
1  m.  From  Fig.  I- 10,  the  optimum  experimental  condition  may  be  determined. 


Fig.  1-10.  Equivalent  radius  of  unbloomed  beam  ^ 
for  six  aperture  sizes.  Gaussian  beam  is  truncated  ^ 
at  1/e^  power  level  and  has  wavefront  curvature  I 
of  1-m  radius.  All  curves  apply  to  wavenumber 
k=  1.2X105  cm"''. 


Figure  I- 11  shows  the  5-mm-aperture  beam  with  thermal  blooming  (curve  B).  By  applying 
an  initial  phase  correction,  thermal  blooming  is  considerably  reduced  (curve  C),  approaching 
the  unbloomed  condition  (curve  A). 

Figure  1-12  shows  the  equivalent  curves  for  a  2 -mm -aperture  beam. 

J.  Herrmann 
L.  C.  Bradley 
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Fig.  1-11.  For  a  beam  with  a  5-mm-diameter 
aperture,  curve  A  is  the  unbloomed  case,  curve  B 
is  with  thermal  blooming,  and  curve  C  is  with 
initial  phase  correction.  Blooming  parameters 
are:  P  =  1  W,  a  =  1 0"^  cm“^ ,  v  =  0.2  cm/sec, 
and  (d  =  0.2  rad/sec. 


Fig.  1-12.  Same  as  Fig.  1-1 1  for  2-mm-aperture 
beam.  Blooming  parameters  also  same,  except 
V  =  0.5  cm/sec  and  w  =  0.5  rad/sec. 
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Fig.  1-13.  Breakdown  thresholds  measured  for  spot 
diameters  0.3  and  0.5  mm.  Data  for  smaller  spot 
size  are  from  TEA  laser  results  at  Lincoln  Labora¬ 
tory  and  from  Ref.  10. 


FOCAL  DIAMETER  (cm) 


Section  I 


B.  EFFECTS 

1.  Laser -Induced  Air  Breakdown 

Several  experiments  were  performed  to  investigate  the  effects  of  atmospheric  dust  on  the 
breakdown  threshold  of  10.6-fmn  laser  beams  with  large  focal  spots. 

6  7 

a.  A  ERL  Breakdown  Experiment  ’ 

In  March  1972,  we  conducted  a  series  of  experiments  at  Avco  Everett  Research  Laboratories 
(AERL)  using  their  1-kJ  pulsed  e-beam  laser.  This  laser  emitted  a  10-  x  20-cm  beam  which 
was  from  i  to  2X  diffraction-limited  and  had  a  peak  pulse  power  k4  X  lO^W.  ^  With  25-m  focusing 
optics,  we  were  able  to  obtain  breakdown  at  a  maximum  spot  area  of  «0. 3  cm  .  The  focused 
power  was  varied  by  either  changing  the  sustainer  voltage  or  by  inserting  finely  perforated 
screens  with  geometrical  transmission  of  0. 70  and  0. 5  in  the  beam.  The  threshold  for  breakdown 
was  obtained  to  within  a  factor  of  approximately  two.  The  results  for  large  spot  sizes  are  shown^ 
in  Fig.  1-13  together  with  Lincoln  Laboratory  TEA  laser  results  and  the  results  from  D.  C.  Smith. 
As  can  be  seen,  the  threshold  continues  to  decrease  at  these  spot  sizes. 

In  general,  for  flux  densities  larger  than  the  threshold,  breakdown  was  not  confined  to  the 
focal  volume  but  instead  consisted  of  a  string  of  discrete  breakdowns  extending  up  the  beam 
(Fig.  1-14).  The  length  of  this  string  increased  with  increasing  laser  power,  as  shown  in  Fig.  1-15. 
The  breakdown  string  always  extended  well  beyond  the  point  where  the  beam  flux  was  equal  to  the 
clean-air  breakdown  value  (»3  x  10^  W/cm^).  The  minimum  flux  densities  where  breakdown  oc- 
7/2 

cur  red  were  W/cm  . 

These  results  can  be  understood  by  assuming  the  discrete  breakdowns  are  occurring  off 
large  dust  particles  in  the  beam.  The  size  distribution  of  dust  in  the  experimental  area  was 
monitored  using  a  light  scattering  particle  counter  (Royco  Model  225)  and  a  256-channel  pulse- 
height  analyzer  (Northern  Model  NS-633).  Particles  in  the  size  range  1  to  12  pm  could  be  meas¬ 
ured  in  this  way.  A  typical  measured  distribution  is  shown  in  Fig.  I-l6. 


15 


dn/dD  (cm~^) 


Section  i 


Fig.  1-15.  Maximum  extent  of  breakdown  string  as  a  function  of  peak  laser  power. 


Fig.  1-16.  Measured  dust-particle-size 
distribution. 
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Fig.  1-17.  Beam  geometry  and  notation.  Breakdown  string  extends  to  position  Z 
where  flux  density  is  l(Z^). 

The  functional  dependence  of  the  breakdown  threshold  on  particle  size  can  be  obtained  in  the 
following  way:  Let  ihe  flux  density  at  the  furthest  extent  of  the  breakdown  string 

(Fig.  1-17).  Assume  that  the  breakdown  is  determined  by  the  largest  particle  size  that  can  be 
present  near  this  position.  The  density  of  particles  of  size  or  larger  is  given  by 

N{D  )  =  C  ^  dD  10'^^  . 

O  Jp) 

o 

The  largest  particle  that  can  be  present  near  is  given  by  the  conditions 
N(D^)  V(Z^)^10 

where  V(Z  )  is  the  volume  of  the  laser  beam  from  the  focal  spot  to  Z  .  That  is,  there  will  be  a 
o  ^  o  ' 

particle  of  size  near  Z^  if  there  are  approximately  ten  such  particles  in  V(Z^),  This  gives 
which,  together  with  1(2-^),  gives  the  threshold  intensity  on  particle  size  shown  in 
Fig.  1-18.  Note  that  the  results  for  are  not  very  sensitive  to  the  number  of  particles  assumed 
to  be  present  in  V(Z^). 

We  conclude  that  the  experimental  results  for  breakdown  threshold  can  be  understood  if  the 
breakdown  threshold  depends  on  aerosol  particle  size  in  the  manner  shown  in  Fig.  1-18. 


Fig.  1-18.  Particle-size  dependence  of  breakdown 
threshold  inferred  from  measured  values  of  zm 
and  dn/dD. 
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Breakdown  Results  with  TEA  Lasers 


It  has  been  reported  recently  that  the  breakdown  threshold  for  10,6~fxm  radiation  decreases 

with  increasing  spot  size  well  below  the  value  predicted  by  classical  microwave  breakdown 

-  2  -1 

theory.  We  have  verified  the  effect  for  focal  spot  sizes  in  the  range  10  to  1  0  cm,  but  found 
that  the  threshold  was  sensitive  to  the  amount  of  dust  near  the  focal  volume,  especially  at  large 
spot  sizes. 

A  chamber  was  designed  and  constructed  which  provided  a  variable  atmosphere  from  dusty 
to  highly  filtered  clean  air.  The  dust  count  and  size  distribution  within  the  chamber  were  mon¬ 
itored  with  a  light  scattering  particle  counter  (Royco  225)  and  a  256-channel  pulse-height  analyzer 
(Northern  NS-633).  We  were  thus  able  to  measure  the  size  distribution  in  the  range  1  to  12  jim. 

A  series  of  threshold  measurements  was  conducted  in  known  atmospheres.  For  the  very- 
clean-air  case,  the  threshold  approached  an  as3rmptotic  value  at  large  spot  size  which  was  close 
to  that  expected  from  microwave  breakdown  theory.  Figure  1-19  shows  these  data  along  with  a 
theoretical  curve  obtained  from  Ref.  9.  The  focused  spot  was  rectangular  with  dimensions 
S  X  5S.  The  fit  to  the  data  was  for  a  diffusion  length  A  =  S/6.  The  other  data  in  Fig.  1-19  corre¬ 
spond  to  an  atmosphere  of  carbon  dust.  As  can  be  seen,  the  threshold  continues  to  decrease 
with  increasing  spot  size.  The  relative  dust  conditions  for  the  three  cases  are  indicated  in  the 
figure  by  the  total  density  of  particles  with  diameters  larger  than  2  jam. 

As  in  the  previous  section,  these  results  can  be  explained  by  assuming  that  the  breakdown 
threshold  is  determined  by  the  largest  particle  present  in  the  focal  volume.  Thus,  for  example, 
at  S  =  0.03  cm  in  Fig.  1-20  the  three  dust  conditions  give  maximum  particle  sizes  of  3,  1,  and 
«1  fxm  for  the  three  data  points.  The  inferred  threshold  dependence  on  particle  size  obtained  in 

this  way  is  shown  in  region  A  of  the  figure.  Also  shown  are  the  results  discussed  in  the  previous 

1 0 

section  and  the  thresholds  measured  directly  for  thin  glass  fibers  and  single  particles. 

The  above  results  indicate  that  atmospheric  aerosols  can  account  for  the  threshold  lowering 
at  large  spot  sizes.  Preliminary  calculations  of  the  cascade  process  in  exploding  particle  vapors 
give  the  right  order  of  magnitude  as  well  as  the  particle- size  dependence  of  the  breakdown 
threshold. 

Further  work  in  this  area  will  include  studies  of  different  material  aerosols,  some  refine¬ 
ment  of  the  theoretical  models,  and  studies  of  the  interaction  of  radiation  with  aerosols  below 


the  breakdown  threshold. 


D.  E.  Lencioni 
L.  C.  Marquet 


2.  Surface  Interactions 

a.  Impulse  Delivery  by  a  Pulsed  Laser  at  10.6  (am 

Measurements  and  a  correlation  of  impulse  delivered  to  targets  of  various  sizes  and  ma¬ 
terials  with  the  use  of  the  A  ERL  1-kJ  CO^  laser  are  reported  in  this  section.  Some  of  these 

11 

measurements  and  the  experimental  procedure  are  reported  in  detail  elsewhere. 

The  total  impulse  delivered  to  carbon  and  aluminum  targets  of  different  sizes  was  measured 

2  12 

for  focal  spot  areas  of  approximately  0.1  and  0.3  cm  (27-m  focal-length  mirror).  AERL  data 

2 

for  aluminum  and  tungsten  targets  are  also  reported  for  focal  spot  areas  of  0.02  and  0.6  cm 
(10.6-m  focal-length  mirror).  The  delivered  impulse  was  found  to  be  independent  of  the  target 
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material  (carbon,  aluminum,  and  tungsten).  There  was  no  measurable  mass  loss  from  the 

aluminum  targets  (Am  <  lO'^^g)  and  the  mass  loss  from  the  carbon  targets  was  less  than  10  g. 

2 

The  projected  target  areas  ranged  from  0.7  to  40  cm  . 

A  simple  model  based  on  an  expanding  cylindrical  blast  wave  and  a  plasma  front  velocity 
proportional  to  the  cube  root  of  the  power  density  is  used  to  obtain  impulse  scaling  relations.  It 
is  shown  in  Refs.  13  and  14  that  the  cylindrical-blast-wave  relations  can  be  modified  to  approx¬ 
imately  account  for  time-dependent  energy  deposition.  The  shock-wave  radius  and  the  resulting 
pressure  are  shown  to  be  proportional  to  the  following  relations: 


R(t)  oc 


(1-5) 


P(t)  oc  E^(t) 


^  [E,(t)]^/^  dt 


-1 


(1-6) 


where  the  energy-per-unit-length  is  defined  as 


E(t) 
I  -  i(t) 


(1-7) 


Here,  E(t)  is  the  total  energy  deposited  in  the  plasma  and  i(t)  is  the  distance  the  plasma  front  has 
moved  away  from  the  target  during  the  time  t.  Then,  for  the  case  of  a  constant  power  density, 
we  find 


2/3 


(1-8) 


where  is  the  focal  spot  area  and  q  is  the  average  power  density  incident  on  the  target.  The 
impulse  is  given  by 


I  <x  P(t)  R(t)^  dt  (1-9) 

and  for  the  case  of  a  constant  power  density 

I  oc  A  ( T  .  (I-IO) 

o  ^ 

Two  limiting  cases  may  be  identified  which  determine  the  termination  of  the  impulse  delivery  r. 
For  small  free-standing  targets,  the  impulse  delivery  ends  when  the  cprlindrical  shock  wave  ex¬ 
pands  off  the  target.  In  this  case,  Eq.  (1-5)  yields  r  ^  ^  (q  )  ,  where  A^  is  the  target 

area.  The  impulse  scaling  relation  is 

ICC  A^Ay2(q)l/3  _ 

For  an  infinite  target,  the  impulse  delivery  ends  when  the  shock- wave  pressure  has  decreased 
to  about  1  atm.  Then,  Eq.  (1-6)  yields  r  A^/^(q)^/^  and  the  impulse  scaling  relation  becomes 

IccA^3/2(_)^^y2^^J  (1.12) 

where  E  is  the  total  energy  in  the  laser  pulse  and  is  the  pulse  length. 
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The  impulse  normalized  by  the  target  area  and  the  square  root  of  the  focal  spot  area  are 
presented  in  Fig.  1-21  as  a  function  of  the  average  power  density  as  suggested  by  Eq.  (I- 11 ).  The 
solid  line  is  a  least-squares  fit  to  the  data  and  has  a  slope  of  0.27.  The  dashed  line  is  fit  to  the 
data  with  a  slope  of  l/3  as  predicted  by  Eq.  (I-ll ).  The  scatter  of  the  data  makes  it  impossible 
to  determine  the  precise  power  of  the  impulse  dependence  on  the  average  power  density,  but  a 
small  exponent  is  indicated  (i.e.,  l/5  to  l/3).  It  should  be  noted  that,  as  pointed  out  in  Ref.  13, 
the  power  density  is  only  known  to  within  about  a  factor  of  3.  Equation  (I-ll )  indicates  that  the 
normalized  impulse  is  relatively  insensitive  to  the  power  density  and  correlates  the  data  to 
within  about  a  factor  of  2. 


Fig.  1-21 ,  Plot  of  normalized  impulse  as  a  function  of  overage  power  density  for  doto  obtained 
with  use  of  AERL  e-beam  CO2  laser.  Solid  symbols  denote  AERL  data  from  Ref.  14.  Solid  line 
is  least-squares  fit  to  data,  1/A^  =  2(q  )0.27.  dashed  line  is  fit  to  data  with  slope  of  1/3. 


Increasing  the  total  energy  above  300  J  led  to  the  formation  of  multiple  aerosol-induced 

breakdowns  in  the  laser  beam.  An  open- shutter  photograph  of  this  effect  is  shown  in  Fig.  1-22. 

2 

The  total  impulse  delivered  to  an  effectively  infinite  target  (40  cm  )  was  measured  on  a  number 
of  runs  both  with  and  without  the  formation  of  aerosol- induced  breakdowns.  In  Fig.  1-23,  the 
measured  impulse  is  shown  as  a  function  of  total  energy  in  the  laser  pulse.  The  limited  data 
show  a  linear  increase  of  impulse  with  energy  in  agreement  with  Eq.  (1-12)  when  multiple  break¬ 
downs  were  not  observed,  and  an  apparent  saturation  of  delivered  impulse  when  multiple  break¬ 
downs  were  observed.  This  interpretation  of  the  data  indicates  that,  with  AERL  machine  pulse 
shape  (an  initial  large  spike),  the  increased  energy  above  3  00  to  400  J  is  being  absorbed  in  the 
breakdowns  and  is  decoupled  from  the  target.* 


*  An  alternative  interpretation  of  these  data  is  possible  if  the  assumptions  of  impulse  delivery  by  a  spherical  blast  wave 
(instantaneous  deposition  of  energy  at  a  point)  and  a  finite  target  are  made.Hjl^  This  model  predicts  that  the  impulse  is 
proportional  to  the  square  root  of  total  energy.  The  impulse  data  shown  in  Fig.  1-23  could  be  fit  to  within  experimental  error  by  a 
function  of  the  square  root  of  the  energy.  This  interpretation  then  says  the  aerosol-induced  breakdowns  absorb  only  a  small  amount 
of  energy  and  do  not  affect  the  measured  impulse.  It  is  clear  that  more  data  are  necessary  to  really  settle  this  point. 
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b.  Plasma  Properties  Associated  with  Irradiation  of  a  Surface  by  a  Pulsed  Laser 

A  number  of  experiments  were  performed  with  the  Avco  device  to  try  to  quantify  some  of  the 
properties  of  the  plasma  produced  in  the  laser-target  interaction.  The  results  of  some  earlier 
camera  data  (Ref.  11)  indicated  that  a  self-luminous  front  with  velocities  of  2  X  10  to  10  cm/sec 
moved  back  in  the  direction  of  the  laser  beam.  To  supplement  this  information  and  provide  two^ 
dimensional  data,  the  streak  camera  was  replaced  by  a  framing  camera  running  at  about  2X10 
frames/sec.  In  addition,  the  plasma  was  back- illuminated  by  a  40-p.sec-long  pulse  of  4881  radi¬ 
ation,  so  that  essentially  the  second  derivative  of  index  of  refraction  (shadowgraph)  rather  than 
the  self-lumination  was  viewed.  The  target  was  a  carbon  surface  and  the  incident  laser  energy 
was  approximately  200  J,  Some  typical  results  are  shown  in  Figs.  1-24  through  1-26. 

Figure  1-24  shows  the  growth  of  the  plasma  on  the  target.  The  growth  is  approximately  that 
of  a  spherical  wave,  with  some  asymmetry  resulting  because  the  impinging  laser  energy  tended 
to  arrive  in  two  separate  cones  rather  than  in  one.  Figure  1-25  shows  similar  data,  except 
that  in  this  case  a  fine  copper  wire  was  placed  in  front  of  the  target,  and  breakdown  initiated  on 
the  wire  shows  a  distinctly  spherical  expansion  with  an  expansion  velocity  of  1.8  X  10  cm/sec. 
Finally,  Fig.  1-26  shows  the  resulting  pattern  in  the  late  stages  of  expansion,  after  the  laser  pulse 
is  completed  and  the  initial  shock  wave  has  moved  away.  Notice  the  turbulent  state  of  the  medium. 
By  using  these  types  of  data  and  neglecting  the  asymmetry  due  to  the  laser  intensity  distribution, 
the  lateral  and  radial  growth  of  a  plasma  ball  was  measured  and  is  shown  in  Fig.  1-27. 

These  data  were  used  to  help  determine  the  absorptivity  of  the  plasma.  In  an  experimental 
arrangement  similar  to  that  described  in  Ref.  11  and  illustrated  in  Fig.  1-28,  the  transmission 


0 

Fig.  1-24. 
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Back-illuminaled  framing  camera  data 
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showing  growth  of  a  plasma  on  a  target. 
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Fig.  1-25.  Back-il iuminated  framing  data  showing  initiation  and  growth  of  a  plasma 
on  a  fine  copper  wire.  Interframe  time  is  0.5  psec. 
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Fig,  1-26.  Late  time  development  of  a  plasma  Initiated  on  a  surface.  Note  turbulent  nature 
of  plasma  medium.  These  data  are  typical  '^30  psec  after  initiation  of  10-pm  laser  pulse. 
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TIME  (/tsec) 

Fig.  1-27.  Lateral  and  forward  growth  of  plasma  on  a  carbon  surface.  Energy  in  pulse 
was  about  200  J, 


10.6->tm  LASER 


Fig-  1-28,  Experimental  arrangement  used  to  measure  plasma  opacity.  An  ion-argon  laser 
(X  =  4881  A)  and  an  order  of  10-pm  high-power  laser  are  made  to  traverse  same  path  through 
plasma  medium  a  distance  d  in  front  of  a  target.  M]  and  M2  are  bending  mirrors,  D]  and 
D2  are  Ge  photon-drag  detectors,  D3  is  a  photodiode  detector,  and  F]  is  a  chain  of  narrow- 
band  4881 -A  filters. 
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history  of  10.6-|jim  and  4881-A  radiation  through  the  plasma  was  monitored  a  distance  d  in 
front  of  the  target. 

The  10.6-}Jim  radiation  used  to  probe  the  plasma  was  obtained  by  taking  one  of  the  grating 
side  orders  and  transmitting  it  through  the  plasma  and  into  a  photon-drag  detector  D^.  At  the 
same  time,  the  unattenuated  laser  pulse  was  monitored  with  photon- drag  detector  D^.  A  1-W 
40-(jisec-long  pulsed  argon-ion  laser  was  triggered  approximately  lOpsec  before  the  CO^  laser 
pulse  and  passed  through  the  same  path  as  the  deflected  side  order  of  the  CO^  pulse.  The  argon 
laser  pulse  was  passed  through  a  series  of  4881-A  narrowband  filters  into  a  Si  photodiode  D^. 


Fig.  1-29.  Normalized  10-pm  pulse  shapes 
a  different  distance  in  front  of  a  carbon  tar¬ 
get.  Solid  curves  are  unattenuated  pulses, 
while  dotted  curves  represent  pulses  after 
being  transmitted  through  plasma  medium, 
(a),  (b),  and  (c)  represent  data  taken  at  dis¬ 
tances  of  0,  3,  and  10  mm  in  front  of  tar¬ 
get,  respectively. 
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Measurements  were  taken  at  three  distances:  0,  3,  and  10mm  in  front  of  the  target.  For 
the  10-pm  data,  the  two  pulses  (attenuated  and  unattenuated)  were  normalized  to  peak  unit  inten¬ 
sity.  The  results  are  shown  in  Fig,  I-29(a-c),  Notice  that  the  CO^  laser  pulse  has  a  rise  time 
of  400  to  500 nsec.  Absorption  begins  after  a  delay  which  is  consistent  with  velocities  in  the 
range  3  x  10^  to  10^  cm/sec.  The  framing  camera  data  were  used  to  determine  the  optical  path 
length  5(d,  t)  at  a  given  distance  d  in  front  of  the  target  and  at  a  given  time.  The  absorption 
coefficient 

Inl/I 

o 


was  then  calculated  and  is  shown  in  Fig.  1-30.  The  simplifying  assumption  of  a  uniform  absorp¬ 
tion  across  the  plasma  was  used. 
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Data  from  the  4881-A  laser  are  shown  in  Figs.  I-31(a-d).  The  relative  transmission  at  d  =  0, 

3,  and  10  mm  are  shown  in  Figs.  I-31(a-c),  respectively.  Figure  I- 31(d)  is  taken  at  0mm,  but 
one  of  the  narrowband  filters  in  front  of  detector  was  removed  so  that  the  plasma  self-emission 
is  shown.  Notice  that  there  is  a  narrow  absorption  spike  (<1  |xsec)  before  the  self-emission  of  the 


plasma  builds  up. 

By  analyzing  these  data  in  the  same  way  as  the  lO-jJim  data,  the  absorption  coefficient  at 
4881  A  is  shown  in  Fig.  1-32, 


Fig.  1-32,  Time-dependent  absorption  coefficient 
data  at  4881  A  at  three  distances  in  front  of  target. 
Because  of  jitter  in  firing  of  this  laser,  uncertainty 
in  initial  points  of  these  curves  is  psec. 


If  one  makes  the  assumption  that,  for  the  10-pm  radiation,  the  absorption  is  due  to  inverse 
bremsstrahlung,  one  can  use  the  absorption  coefficient  to  determine  the  electron  density,^ 
using  the  formula 


o!(cm"^)  =  1.6  X  10  -  e 


-hv/kT  kT^ 


') 


hr 


(1-13) 


Here,  a  is  the  absorption  coefficient  in  cm~^,  \  is  the  wavelength  in  pm,  r  is  the  electron 
density  in  number /cm  ,  is  the  electron  temperature  in  eV,  and  r  is  the  laser  frequency. 

The  electron-density  histories  shown  in  Fig.  1-33  were  calculated  with  Eq,  (1-13),  the  absorption 
coefficients  given  in  Fig.  1-30,  and  the  assumption  of  a  1-eV  plasma.  Also  shown  in  Fig.  1-33 
are  electron  densities  obtained  in  the  same  manner  from  the  4881-A  absorption  data.  The  as¬ 
sumption  of  a  1-eV  plasma  temperature  is  a  reasonable  estimate  of  the  actual  temperature; 
thus,  the  calculated  electron  densities  are  of  the  right  order  of  magnitude  and  give  correct  rel¬ 
ative  values. 
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Fig,  1-33.  Electron  densities  calculated  from  data  of  Figs.  1-30  and  1-32  for  1 -eV  temperature 
and  assumption  that  predominant  absorption  mechanism  is  inverse  bremsstrah lung.  This  assumption 
is  plausible  for  10-pm  data  and  gives  reasonable  estimate  for  electron  density  time  history  in 
1017  to  3  X  10l7  cm-3.  The  4881 -A  data  suggest  some  photo-excited  absorbing  species  as  more 
probable  candidates  for  vapor  opacity. 
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These  data  indicate  that  the  absorption  at  4881  A  is  not  primarily  due  to  free  electrons.  In 

O 

fact,  to  within  about  ~1  |jLsec,  the  absorption  at  4881  A  begins  simultaneously  at  all  three  dis¬ 
tances  d  in  front  of  the  target.  Recall  the  absorption  dip  at  the  beginning  of  the  10-pm  laser 
pulse  in  Fig.  I- 31(d)  and  notice  that  the  absorption  takes  place  before  there  is  appreciable  self¬ 
emission  from  the  plasma,  i.e.,  the  plasma  is  still  relatively  cool  with  no  appreciable  electron 
density. 

The  simultaneous  nature  of  the  absorption  onset  suggests  the  production  of  a  photo-excited 
species  which  absorbs  strongly  at  4881  A  and  decays  with  a  lifetime  of  about  lOpsec. 

H.  Kleiman 
R.  W.  O'Neil 

3.  Explosive  Vaporization  of  Fog  Droplets  by  a  10.6-juim  Laser  Pulse 

Study  continues  of  the  explosive  vaporization  of  fog  droplets  by  a  10.6-fmi  laser  pulse.  The 
early  stage  (up  to  about  2  psec)  of  the  explosion  process  was  described  in  previous  Optics  Re¬ 
search  Reports  *  We  studied  the  later  stage  of  the  process  by  observing  the  vapor  "blob"  pro¬ 
duced  by  the  evaporation  of  single  droplets.  These  observations  were  made  with  a  high-speed 
schlieren  photography  system  that  was  described  in  the  last  Optics  Research  Report  [(1971:2), 
DDC  AD-901213];  early  results  on  the  vaporization  of  oversize  droplets  were  also  presented  in 
that  report. 

In  a  ground  fog  or  a  cloud,  the  water  droplets  usually  range  in  size  from  <1-  to  ~25-p-m  ra¬ 
dius.  In  the  present  study,  droplets  ranging  from  5-  to  25-pjn  radius  have  been  irradiated.  In 
addition,  some  oversize  droplets  (25-  to  70-pm  radius)  have  been  studied  because  their  behavior 
does  help  in  understanding  the  processes  involved.  The  laser  energy  at  the  position  of  the  water 
droplet  was  varied  between  5  and  15  j/cm  in  these  experiments.  The  peak  power  density  was 
typically  10  to  30MW/cm  . 

Electrical  breakdown  did  not  occur  under  these  conditions.  This  was  determined  by  means 
of  time-exposed  photomicrographs  of  the  exploding  droplet;  these  photos  were  taken  without  any 
back- illumination.  No  light  was  emitted  during  the  explosive  vaporization  of  a  water  droplet. 

If,  however,  instead  of  a  water  droplet,  a  25-fJLm-radius  glass  fiber  is  placed  at  the  same  posi¬ 
tion  and  irradiated  with  an  identical  pulse,  electrical  breakdown  does  occur  and  a  brilliant  flash 
of  light  is  produced. 

a.  Results  and  Discussion 

The  schlieren  pictures  show  an  expanding  sphere  of  vapor  surrounded  by  a  rapidly  expanding 
shock  wave.  Actually,  the  observed  "vapor  sphere"  is  a  sphere  of  air  that  has  been  heated  by 
the  vapor.  The  sphere  of  heated  air  has  an  index  of  refraction  different  from  that  of  the  ambient 
medium,  and  this  difference  is  sensed  by  the  schlieren  system. 

Figure  1-34  shows  schlieren  photographs  of  four  exploding  droplets  at  times  from  0.7  to 
37  fjLsec  after  the  10.6-fjLm  laser  pulse.  The  radius  of  each  water  droplet  is  indicated  in  the  pic¬ 
ture.  At  0.7-}jLsec,  the  vapor  sphere  cannot  be  distinguished  within  the  shock  wave,  but  the 
second  photograph  (4.9  psec)  clearly  shows  tlie  vapor  sphere  at  the  center  of  the  shock-wave 


*  Optics  Research  Reports,  Lincoln  Laboratory,  M.LT.  (1970:3),  DDC  AD-882617,  and  (1971:1),  DDC  AD-888823-1 
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|_-  5  -  3  612-21 


25  pm  16. 2Msec  31  pm  37ps€c 


Fig.  1-34.  Schlieren  photographs  of  four  exploding  water  droplets.  Radius  of  each  droplet 
is  shown  at  left,  and  elapsed  time  since  10.6-pm  laser  pulse  is  shown  at  right.  Laser  pulse 
direction:  left  to  right. 

sphere.  In  the  third  picture  (16.2  psec),  the  shock  wave  can  be  seen  at  the  top,  after  it  has  re¬ 
flected  from  the  horizontal  surface  (which  happens  to  be  the  bottom  of  the  droplet  generator). 
Finally  (at  37psec),  the  shock  wave  is  beyond  the  FOV  of  the  system.  This  sequence  of  pictures 
also  shows  the  expanding  volume  of  the  vapor  sphere.  In  these  photographs,  the  direction  of  the 
10.6-pm  laser  pulse  is  from  left  to  right;  the  schlieren  knife-edge  is  oriented  horizontally.  (The 
pictures  also  show  some  extraneous  features  that  are  common  to  all  the  photographs;  these  are 
due  to  aberrations  in  the  optical  system  and  to  imperfections  in  the  lenses.  They  can  be  readily 
distinguished  and  disregarded. ) 

In  the  study  of  the  early  stage  of  the  water-droplet  explosion,  we  showed  that  the  liquid-vapor 
mixture  was  expanding  at  a  radial  velocity  of  0. 13  .mm/psec  in  the  time  interval  0  to  2.3  psec. 

Our  schlieren  photographs  show  that  the  vapor  sphere  expands  with  a  radial  velocity  of  0.015  mm/ 
psec  in  the  time  interval  5  to  16  psec,  and  the  velocity  has  decreased  to  a  value  of  0.  OlOmm/psec 
in  the  interval  16  to  37  psec.  Even  in  this  interval,  the  radial  expansion  rate  is  ~300  times 
greater  than  the  diffusion  velocity,  indicating  that  the  expansion  process  is  still  dominated  by 
the  momentum  of  the  initial  explosion. 

In  Fig.  1-35,  the  photographs  show  the  effect  of  droplet  size  on  the  magnitude  of  the  observ¬ 
able  effects.  The  small  droplets  (11-  and  13-pm  radius)  at  the  left  produced  a  faint  but  clearly 
observable  shock  wave  with  a  faintly  distinguishable  vapor  sphere  at  the  center.  The  pictures 
at  the  right  show  the  much  stronger  effects  produced  by  the  explosive  vaporization  of  25 -pm- 
radius  droplets,  where  the  mass  of  water  is  ~8  times  as  great.  In  each  of  these  shots,  the  energy 
density  of  the  laser  pulse  was  about  10  j/cm  ,  and  the  pulse  shape  was  about  the  same. 


32 


Section  I 


Fig.  1-35.  Schlieren  photographs  of  four  exploding  water  droplets,  showing  effect  of  droplet 
size.  Laser  pulse  direction:  left  to  right. 


Under  certain  conditions,  the  hot  vapor  assumed  a  nonspherical  shape  —  usually  when  the 
water  droplet  was  very  large  (>  35-pm  radius)  or  when  the  energy  density  was  very  low  (~5  J/ cm  ). 
In  these  cases,  the  vapor  volume  was  shaped  like  a  dumbbell  with  its  axis  oriented  approximately 
along  the  direction  of  the  laser  pulse.  Figure  1-36  shows  a  typical  example,  an  ~35-|JLm-radius 
droplet  at  6.7psec  after  the  laser  pulse.  The  left  sphere  of  the  dumbbell  is  at  the  center  of  the 
shock  wave,  and  this  is  the  site  of  the  original  explosion.  If  the  absorbed  energy  is  not  sufficient 
to  cause  immediate  (within  psec)  vaporization  of  the  droplet,  and  if  the  droplet  is  large  enough 
to  undergo  surface  heating  rather  than  volume  heating,  then  the  remainder  of  the  evaporating 


Fig.  1-36.  Schlieren  photograph  of  large 
(>35-pm  radius)  exploding  water  droplet 
6.7  psec  after  10.6-pm  laser  pulse.  Laser 
pulse  direction:  left  to  right. 
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droplet  is  propelled  by  the  rocket  action  of  its  vaporizing  surface.  The  motion  is  always  oppo¬ 
site  to  the  direction  that  the  blow-off  vapor  is  moving,  as  required  to  conserve  momentum. 

A  picture  similar  to  Fig.  1-36  was  obtained  with  a  26 -jam- radius  droplet  irradiated  with  a 
laser  pulse  of  4  J/cm^.  But  when  another  droplet  of  the  same  size  was  irradiated  with  a  9-J/cm 
pulse,  the  usual  spherical  vapor  volume  resulted.  This  behavior  can  be  explained  by  the  curves 
in  Fig.  1-37,  which  show  the  energy  absorbed  by  a  water  droplet  for  pulses  of  various  energy 
densities.  Actually,  this  is  the  energy  that  would  be  absorbed  by  a  water  droplet  that  remained 
intact  during  the  laser  pulse.  These  curves  were  drawn  using  the  absorption  cross  sections  in 
the  tables  of  Zel'manovich  and  Shifrin.^'^  In  the  case  of  the  9- J/cm  pulse,  the  energy  that 


Fig.  1-37.  Energy  absorbed  by  wafer  droplet  for  incident  10.6-pm  laser  pulses  of  indicated 
energy  density.  Energy  of  vaporization,  is  also  shown  for  comparison. 

could  be  absorbed  by  the  water  droplet  is  substantially  greater  than  the  energy  required  for 
vaporization.  In  the  case  of  the  4- j/cm^  pulse,  the  energy  absorbed  by  the  droplet  is  much  less 
than  the  energy  of  vaporization.  Our  early  studies  using  the  tri-scopic  camera  system  showed 
that  the  droplets  of  radius  >  15  pm  undergo  surface  heating.  Such  droplets  when  partially  vapor¬ 
ized  can  be  propelled  along  the  general  direction  of  the  laser  pulse;  this  is  the  case  for  the  drop¬ 
let  in  Fig.  1-36,  which  could  absorb  barely  enough  energy  for  vaporization.  In  general,  the 
asymmetric  heating  (surface  heating)  of  optically  thick  droplets  will  result  in  an  asymmetric 
vapor  volume  unless  the  absorbed  energy  is  much  greater  than  that  required  for  vaporizing  the 
droplet.  The  asymmetric  volume  will  be  shaped  like  a  dumbbell  or  like  an  incipient  dumbbell, 
with  the  second  sphere  undeveloped  or  partially  developed. 
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b.  Energy  of  the  Explosion 

Figure  1-38  shows  the  radius  of  the  shock  wave  as  a  function  of  time  for  many  different 
shots,  covering  three  ranges  of  droplet  size.  Although  the  points  fall  close  to  the  line  of  slope 
corresponding  to  a  shock- wave  velocity  of  Mach  1  (in  air),  it  will  be  noted  that  the  zero- time 
intercept  is  0.5mm,  which  is  20  times  greater  than  the  initial  droplet  radius.  In  the  first  couple 
of  microseconds,  the  shock-wave  velocity  is  much  greater  than  Mach  1. 


Fig.  1-38,  Radius  of  shock  wave  as  function  of  time. 

The  radius  of  the  shock  wave  at  a  given  time  can  be  used  to  calculate  the  explosive  energy, 
as  discussed  in  the  last  Optics  Research  Report  [(1971:2),  DDC  AD-901213].  We  have  calculated 
the  explosive  yield  for  5  shots  involving  water  droplets  of  25 -pm  radius.  The  yield  may  be  com¬ 
pared  with  the  difference  between  the  absorbed  energy,  as  determined  from  Fig.  1-37,  and  the 
energy  needed  to  evaporate  the  droplets.  The  agreement  is  within  about  25  percent.  This  en¬ 
ergy  balance  cannot  be  used  for  large  (  >40-pm-radius)  droplets  because  they  do  not  evaporate 
completely. 

P.  Kafalas 
J.  Herrmann 
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Fig.  1-39.  Schematic  diagram  of  500-J  EDL. 
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C.  DEVICES 

1.  Electron  Beam  (e-Beam) 

a.  500-J,  e-Beam  Excited  Laser 

Figure  1-39  is  a  schematic  diagram  of  the  500-J  EDL  now  under  construction*  All  major 
components  are  on  hand  and  assembly  is  nearly  complete.  In  keeping  with  its  planned  use  as 
a  source  of  high-energy  10.6-|jim  pulses  of  good  optical  quality,  the  operating  sequence  has  been 
automated  as  much  as  possible  to  keep  the  output  reproducible  and  reliable. 

The  electron  gun  is  a  triode  and  can  be  grid  controlled  from  cutoff  to  current  densities  in 
excess  of  lOOmA/cm  through  a  15  X  104cm,  1-mil  titanium  foil  window.  It  is  capable  of  pulses 
of  from  1-  to  lOO-jjisec  duration  with  rise  times  of  0.1  psec. 

The  optical  cavity  has  an  effective  length  of  2.5  m  with  a  maximum  cross  section  of  15  x  20cm. 
The  active  volume  is  22.5  liters. 

To  avoid  the  need  for  large,  high- power  windows,  a  fast  shutter  valve  is  opened  to  pass  the 
outcoupled  optical  beam. 

The  optical  cavity  will  be  evacuated  after  each  shot  and  then  refilled  with  the  laser  gas  mix¬ 
ture.  A  gas  blender  will  allow  the  mixture  to  be  easily  optimized  for  the  pulse  length  and  power 
chosen  by  an  experimenter. 

First  operation  in  a  laser  is  expected  early  in  August,  followed  by  about  60  days  of  tests 
and  diagnostics.  We  anticipate  that  the  device  should  be  operational  in  October. 

A.  J.  Morency 
R.  J.  Carbone 


b.  e-Beam  CO^  Laser  Kinetics 

1 8 

The  CO^-N^-He  rate  equations  of  Douglas -Hamilton  were  programmed  using  the  Runge- 
Kutta  method.  In  Douglas -Hamilton's  approach,  the  ^3  inodes  of  CO^  and  the  vibra¬ 

tional  levels  of  are  assumed  to  equilibrate  internally  so  rapidly  that  each  of  these  modes 
can  be  described  by  a  vibrational  temperature.  The  four  equations  describing  these  temperatures, 

plus  a  fifth  equation  describing  the  gas  temperature,  compose  the  coupled  set  of  differential  equa- 

19 

tions.  With  the  electron  excitation  coefficients  calculated  by  Nighan  and  the  relaxation  rates 
reported  by  Taylor  and  Bitterman,^^  the  computed  gains,  pulse  lengths,  and  pulse  energies  were 
in  reasonable  agreement  with  the  experimental  results  on  the  Avco  e-beam  laser. 

In  this  analysis,  the  laser  flux  is  determined  by  the  constraint  that,  during  the  pulse,  the 
laser  gain  is  equal  to  the  cavity  losses.  Accordingly,  these  equations  do  not  correctly  describe 
the  initial  portion  of  the  pulse  where  the  gain  still  greatly  exceeds  the  losses,  and  the  equations 
do  not  predict  the  experimentally  observed  initial  spike.  To  remedy  this,  a  sixth  equation  de¬ 
scribing  the  growth  of  the  stimulated  emission  was  included.  This  indeed  provided  the  initial 
spike,  but  its  computed  50-nsec  width  is  considerably  smaller  than  the  observed  300  to  500 nsec. 
This  discrepancy  may  be  due  to  the  fact  that  on  the  time  scale  of  the  initial  spike,  the  assump¬ 
tion  of  vibrational  temperatures  may  be  invalid.  We  are  presently  investigating  this  possibility. 

S.  Marcus 
D.  Zanni 

*  See  Optics  Research  Report  (1971:2),  p.  23,  DDC  AD-901213. 
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Fig.  1-40.  Geometry  of  500-J  pulse  unstable  confocal  resonator  cavity  for  optimum  mode 
quality-energy  product  (outcoupling  5  =  0.85). 


Fig.  1-41 .  Output  intensity  normalized  to  saturation  intensity  as  function  of  geometric  outcoupling 
[5  =  ]  _(i/rn^)]  for  uniform  mode  filling,  with  mirror  reflectance  =  0.99,  cavity  length  2.5  m, 
active  length  1.0  m,  and  small  signal  gain  0.032cm“^. 
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c.  Cavity  Design  of  a  Pulsed  500- J  Laser 

As  its  first  priority,  the  EDL  effort  has  a  pulsed  CO^  laser  device  that  is  designed  to  out- 
couple  500  J  of  10.6-|JLm  radiation/''  The  volume  of  the  laser  medium  will  be  22.5  liters  con¬ 
tained  within  a  length  of  1  m.  A  2.5-m  unstable,  confocal,  positive  branch  resonator  (shown  in 
Fig.  1-40)  was  chosen  to  effectively  utilize  this  large  volume /short -length  gain  medium. 

The  cavity  design  was  chosen  so  that  optical  energy  will  not  focus  within  the  device,  and  the 
forward  and  backward  energy  flow  will  very  nearly  fill  the  gain  medium.  The  optimized  reso¬ 
nator  parameters  were  chosen  as  a  tradeoff  between  total  power  output  reduction  while  maintain¬ 
ing  a  single  mode  as  the  output  coupling  was  increased.  Figure  1-41  illustrates  the  effect  of 
power  output  vs  coupling  at  a  gain  of  0.  032  cm  ^  and  99-percent  reflecting  mirrors.  The  Rigrod 
approach  was  used  to  determine  the  power  extraction  as  a  function  of  predicted  gain  and  optimal 
energy  storage,  A  reasonable  tradeoff  between  mode  quality  and  optimum  energy  output  is  at  an 
output  calculated  to  be  60  percent.  Using  a  geometrical  approach,  this  out  coupling  will  utilize 
at  least  52  percent  of  the  active  volume  for  an  inscribed  (within  a  15  X  15  cm  square)  circular 
mode,  65  percent  for  a  15  X  15  cm  square  mode,  and  87  percent  for  a  fitted  rectangular  mode  of 
15  X  20  cm.  Several  cavity  outcouplings  will  be  tested  to  experimentally  verify  the  power  output 
predictions,  and  to  compare  the  experimental  results  with  mode  calculations  that  have  been  re- 
ported^^'^^  and  with  calculations  presently  being  done  at  Lincoln  Laboratory. 

The  criterion  used  to  judge  beam  quality  outcoupled  from  the  laser  is  the  energy  deposited 
in  the  Airy  disks  compared  with  81  percent  of  the  energy  deposited  in  the  first  Airy  disk  for  a 
uniformly  illuminated  aperture.  For  the  6  0-percent  coupling  case,  83  percent  of  the  output  en¬ 
ergy  will  be  deposited  within  the  first  three  Airy  disks,  or  3X  diffraction-limited  by  this  criterion. 

Optimization  of  mode  quality  is  achieved  by  maximizing  the  product  of  the  energy  output  and 
the  energy  fraction  in  the  first  Airy  disk,  as  shown  in  Fig.  1-42.  This  optimum  outcoupling  is  at 
85  percent,  which  corresponds  to  a  reduction  of  about  15  percent  in  energy  extracted  per  unit 
volume  (relative  to  60-percent  outcoupling),  approximately  82-percent  volume  filling  factor, 
and  about  45  percent  of  the  net  energy  within  the  first  Airy  disk. 

The  absolute  power  outcoupled  for  the  cavity  design  has  been  calculated  from  the  rate  equa¬ 
tions  and  energy  balance  within  the  laser  medium.  The  approach  used  followed  that  of  Douglas- 

Hamilton.^^  Predictions  of  medium  heating,  pulse  length,  and  gain  are  reliable,  but  calculated 

2  2 

optical  energy  storage /cross- sectional  area  at  7.4  j/cm  with  2mA/cm  e-beam  pumping  may 
be  too  high.  This  indicates  an  extraction  of  74  J/liter  for  a  1-m  device.  Table  I-l  gives  typical 
computed  values  of  pulse  length,  gain,  and  output  energy  for  different  e-beam  current  intensities. 

R.  J.  Carbone 
H.  Granek 


d.  Mode  Calculations  for  the  Unstable  Resonator 

The  Fredholm  integral  equation  for  the  asymmetric  unstable  resonator  cavity  is  used  for 

the  mode  field  distributions  without  a  gain  medium  in  the  mirror  planes.  The  field  distributions 
22 

are  given  by 


*  Optics  Research  Report  (1971:2),  p.  23,  DDC  AD-901213. 
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Fig.  1-42.  Amount  of  output  intensity  (normalized  to  soturation  intensity)  concentrated  within 
first  Airy  dark  ring  of  far-field  pattern  (incident  plane  wave  approximation  onto  outcoupling 
annulus)  for  conditions  given  in  Fig.  1-41. 


TABLE  1-1 

TYPICAL  COMPUTED  LASER  PARAMETERS  FOR  GIVEN  e-BEAM 
CURRENT  INTENSITIES 

Current  Density 

Peak  Zero- Flux  Gain 

Pulse  Length 

Energy  Out 

(mA/cm^) 

(cm"') 

(nsec) 

(J/cm2) 

2 

0. 0282 

43.8 

7.4 

10 

0.0326 

25.7 

10.0 

20 

0.0339 

20.1 

11.1 

30 

0.0345 

17.4 

11.8 
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7^(p^(x)  =  /■'■^(k/L)  yJ^(kxy/L)  expH(k/2L)  (g^x^  +  g^y^)]  ^2^  dy 

y2<?’2<y)  = y  *  zJ^  (kyz/L)  exp[-j(k/2L)  (g2y^  +  g^z^)]  <;d^(z)  dz  .  (1-15) 

In  this  expression,  g.  =  1  -  (L/R^),  where  is  the  radius  of  curvature,  L  is  the  mirror  separa¬ 
tion,  y's  are  eigenvalues  to  be  determined,  and  2a.  is  the  mirror  diameter. 

One  can  substitute  Eq.  (1-15)  into  Eq.  (1-14)  to  determine  an  eigenvalue  equation  for  the  field 

26 

at  the  output  mirror.  Use  of  Bessel  function  identities  permits  elimination  of  the  integral  of 
the  product  Bessel  functions  and  reduction  to  a  form 


dz 


Manipulation  of  this  expression  and  expansion  of  the  unknown  functions  in  the  form 

OO 

<p{x)  exp  [-j  ^  (Zg^g^  -  1 )  x^j  =  >  (x^) 
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(1-16) 


(1-17) 


m=0 


where  =  27ra^^/L  allows  simplification  to  an  infinite  matrix  of  the  form 
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//I  -i_  +  \  I  -i+l  nC./4gj 

{y^y^)  3  ^  ^  exp{-U’  [2j{2g^g^  -  1)  +  1]}  (U> 


IH 


Y  a  L^^*(U') 
^  m  m  '  ^ 


m=0 


(1-18) 


Calculations  of  a  truncated  form  of  this  matrix  for  the  empty  cavity  modes  of  the  unstable  res¬ 
onator  are  under  way.  Comparison  of  these  results  with  propagation  calculations  and  experi¬ 
mental  data  is  anticipated.  Granek 


e.  Longitudinal  e-Beam  Excited  Laser 

The  3-nsec,  600-keV,  e-beam  source  described  in  the  previous  Optics  Research  Report 
was  used  to  ionize  a  CO^-N^-He  medium  from  which  laser  pulses  lasting  up  to  120|jisec  were  ob¬ 
tained.  In  the  configuration  being  used,  the  e-beam  is  nearly  colinear  with  the  laser  axis  and 
normal  to  the  sustainer  electric  field.  An  article  describing  the  operation  of  this  laser  has  been 
published.^^  In  the  meantime,  the  maximum  output  energy  was  increased  to  22.4  J  over  the 
11.5  J  reported  in  that  paper.  This  was  accomplished  by  improving  the  design  of  the  electrodes 
and  by  increasing  the  active  laser  volume.  We  are  presently  redesigning  the  laser  cavity  to 
permit  passive  Q- switching. 

S.  Marcus 
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f.  Pulse  Shaping  a  CO^  Laser  Beam 

For  a  number  of  applications  it  is  desirable  to  be  able  to  shape  in  time  high-energy  laser 
pulses  of  10.6-pm  radiation.  One  approach  is  to  try  to  build  pulse-shaping  capabilities  into  the 
laser  itself,  for  instance,  by  modulating  the  grid  of  an  e-beam  pumped  laser.  However,  for 
other  classes  of  pulsed  machines  it  is  more  practical  to  try  to  develop  external  pulse-shaping 
capability.  Here,  we  report  some  successful  preliminary  attempts  to  generate  simple  waveforms 

using  saturable  absorbers  external  to  the  laser  cavity.  A  similar  technique  has  recently  been 
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reported  for  use  in  the  visible  range. 

The  basic  idea  can  be  understood  by  considering  the  expression  for  absorption  in  a  two-level 
system; 

ri  +  f  exp{-[l  +  (l/lg)]^/^}] 


Here,  ot  is  the  unsaturated  absorption  coefficient,  -  hv/Zar  is  the  saturation  intensity,  o*  is 
the  cross  section  for  molecular  absorption,  v  is  the  laser  frequency,  and  r  is  the  relaxation 
time  of  the  upper  level.  In  all  useful  absorbers,  the  initial  and  final  levels  interact  strongly  with 
closely  lying  rotational  and  vibrational  levels  and,  consequently,  the  two-level  picture  is  not  a 
valid  one.  However,  it  can  be  shown  that,  when  the  relaxation  time  among  these  levels  is  short 
compared  with  r  and  the  laser  pulse  length,  the  above  relation  still  holds  except  that  =  Yiv/2cr6, 
where  the  factor  i/d  is  the  effective  number  of  interacting  levels. 

This  gives  the  medium  saturation  time  constant 


T  = 

s 


s 


Suppose  we  require  the  pulse  length  =  r^.  This  would  allow  the  medium  to  saturate 
gradually  and  generate  a  ramp  function,  which  would  require  an  intensity 


\  _  hK  /  T  _ 

1  ■  2aT6  It  J 

\  p 

/  '  P  ' 

If  the  peak  of  the  ramp  is  to  be  substantially  unattenuated,  then  we  must  have  I  »  or. 


equivalently,  t  t^.  hi  this  case,  a  square  pulse  of  intensity 


I  =  I 


hv 


s  T 


2a  T  6 
P 


should  generate  a  ramp  function.  Figure  1-43  shows  the  experimental  arrangement  in  which  we 
used  a  CO^  laser  which  generated  a  pulse  with  uniform  power  output  of  «4  X  10  W  for  SOpsec 
(see  Ref.  28).  This  pulse  was  focused  through  a  20-cm  absorption  cell  with  Brewster  angle  win¬ 
dows.  A  KCl  flat  diverted  about  5  percent  of  the  energy  to  a  BaTiO^  calorimeter  so  that  the 
pulse  energy  could  be  monitored.  The  pulse  at  the  cell  exit  was  monitored  with  a  photon-drag 
Ge  detector.  The  beam  was  approximately  1  cm^  in  area  at  the  entrance  and  exit  of  the  cell,  and 
perhaps  1  /4  that  area  at  the  center  of  the  cell. 
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Fig.  I“43.  Experimental  arrangement  for  observing  pulse  shaping. 

Figure  I*“44(a-d)  shows  a  series  of  waveforms  generated  at  different  pressures  of  SF^.  In 
the  region  where  a  is  a  function  of  pressure,  the  pressure  and  intensity  can  be  adjusted  to  main¬ 
tain  the  product  la  constant.  This  allows  generation  of  the  same  waveform  at  varying  intensity 
levels. 

Figure  I-45(a-d)  shows  the  results  of  using  Freon  12  as  an  absorbing  gas.  The  relaxation 
time  here  is  short  compared  with  the  pulse  length,  so  that  the  intense  leading  edge  of  the  pulse 
goes  through  unattenuated  while  the  rear  of  the  pulse  is  absorbed.  This  is  the  more  familiar 
case  of  pulse  shaping. 

We  have  demonstrated  the  ability  to  produce  slow  or  fast  opening  attenuators  for  uniform 
laser  pulses,  and  to  use  these  attenuators  to  produce  ramp  and  spikes  of  desired  temporal 

characteristics.  H.  Kleiman 

S.  Marcus 


g.  Feasibility  Study  for  Amplification  of  Short  Pulses 
and  Frequency- Modulated  Pulses 

Here,  we  report  the  results  of  a  feasibility  study  on  the  design  of  a  very  broad-band,  high- 

power  laser  amplification  system.  There  is  current  research  interest  in  an  amplifier  at  10.6  pm 

that  could  handle  either  a  repeated  burst  sequence  of  nanosecond  pulses  or  a  repeated  sequence 

of  pulses  less  than  a  microsecond  in  length  with  FM  up  to  1.5  GHz.  Energy  outputs  in  excess  of 

3  J  per  pulse  and  average  power  of  20  kW  are  required,  with  total  coherence  time  of  the  pulse 

sequence  extending  from  circa  SOpsec  to  5  msec.  To  handle  the  information  bandwidth  of  the 

pulses,  the  amplifier  would  have  to  operate  at  pressures  of  at  least  15  0  to  200  torr  for  5  00- MHz 

30-32 

signals,  and  approaching  an  atmosphere  of  pressure  for  full  bandwidth  capacity.  There  are 

four  devices  in  the  ''high" -pressure  regime  that  could  possibly  satisfy  the  requirements  of  min- 
imum  optical  distortion  (<A./40):  longitudinal  discharge,  cross -flow  laser  ;  waveguide  laser  ; 
preionized  devices;  and  independently  ionized  devices.  The  first  type  has  been  limited  to  pres¬ 
sures  below  100  torr,  thus  does  not  exhibit  sufficient  bandwidth  to  be  applicable.  The  second  is 
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FlgJ-44.  Pulse  shaping  using  SF^  in  absorption  cell.  Pressure  of  SF^  is  (a)  0,  (b)  10, 
(c)20,  and  (d)  35  torr,  respectively. 
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(c)  (d) 


Fig.  1-45.  Pulse  shaping  using  Freon  12  in  absorption  cell.  Pressure  of  Freon-12  is  (a)  0, 
(b)  40,  (c)  100,  and  (d)  200  torr,  respectively. 
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too  small  for  appreciable  energy  extraction  (although  potentially  useful  in  preamplification).  In 
the  third  category,  three  electrode  devices^^’^^  are  a  possibility,  but  the  discharge  optically 
distorts  the  medium  (especially  in  CO,-rich  mixtures)  and  the  gains  are  appreciably  lower  than 
in  the  last  category.  However,  the  e-beam  preionized  laser  '  (see  also  Sec.  C-l-a)  is  a  de¬ 
vice  that  already  approaches  the  capabilities  necessary  for  the  application.  Likewise,  the  UV 
photoionized  laser  presently  being  researched  ^  (see  also  Sec.  C-2-a)  offers  similar  promising 
possibilities  for  this  application. 

Both  in  the  pulse  burst  and  FM  pulse  modes,  the  problems  of  medium  homogeneity  and 
preservation  of  gain  (and  pulse  coherence)  during  variable  time  intervals  are  important,  and 
shall  be  treated  subsequently.  Likewise,  temporal  pulse  distortion  in  both  cases  is  important. 
However,  there  are  differences  that  must  be  treated  separately:  the  ability  of  the  molecular 
medium  to  be  reasonably  utilized  in  times  short  compared  with  its  relaxation  rates,  spurious 
FM  effects  due  to  medium  dispersion,  saturation,  and  pulse  sharpening. 

The  following  relaxation  time  calculations^^  were  done  for  a  HeiN^iCO^  mixture  of  3:2:1  at 
a  total  pressure  of  1  atm.  The  results  are  easily  scaleable  in  pressure.  The  inter-rotational 
relaxation  rate  is  ~0. 12 nsec-atm^^;  therefore,  sufficiently  strong  nanosecond  pulses  will  dump 
all  the  rotational  levels  of  the  001  vibrational  state,  while  such  a  pulse  would  be  too  fast  to  dump 

^  3nsec-atm  (CO^)]  or  allow  replenish- 
12  0 nsec-atm).  Since  these  two  mechanisms  contribute  to  the 


any  of  the  higher  vibrational  00^-^  states  of  CO^  [^qq^  ^  j 


ment  by  the  nitrogen  pump  (r- 

efficiency  of  energy  extraction  in  the  saturation  intensity  would  be  reduced  by  at  least  an 

order  of  magnitude  for  individual  pulses.  However,  the  higher  vibrational  states  would  partially 
repump  before  the  next  pulse  comes  through,  and  the  nitrogen  would  weakly  repump  during  a 
sequence  of  such  pulses. 

A  single  arbitrary  pulse  (in  time)  passing  through  a  high  gain  amplifying  laser  medium  will 
distort  in  shape,  because  its  leading  edge  always  sees  an  unsaturated  gain  and  will  grow  expo¬ 
nentially  (until  other  loss  mechanisms  come  into  play),  while  the  back  sections  of  the  pulse  will 
see  a  successively  more  saturated  gain  medium.  This  results  in  pulse  sharpening  at  the  leading 
edge  (and  a  corresponding  frequency  domain  distortion  that  is  important  in  the  FM  mode).  For 
the  nanosecond  pulses,  the  CO^  laser  acts  as  a  two-level  (total  vibrational  state  populations) 
system  and  the  interaction  of  the  population  inversion  and  optical  field  growth  can  be  shown  to 


be  described  by 


r?(x,  t) 
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exp  I 


[Ncrc  p  (?>j^(t')  dt']  +  exp[-Na  dx' ] ' 


7^(x)  exp[-Na  y^jx)  dx] 


and 


^(x,  t)  = 


(p^{t  -  x/c) 


{l  -  exp[-Nc7  ?]q(x')  dx'])  exp  [-Ncrc  f 


(t-x/c) 


(p^iV)  dt'] 


where 


<P  = 


2$ 


+^2 


V  = 


”2 


N 


"2 
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and 

^  =  number  of  photons /volume 

n^  =  particle  density  in  respective  vibrational  states 

a  =  resonance  absorption  cross  section/active  molecule. 

For  a  uniform  axial  gain  profile  77 (x,  0)  ~  Vq  ^  specified  pulse  form  ^(0,  t),  the  output  pulse 
can  then  be  integrably  determined.  The  distortion  of  a  Lorentzian  pulse  of  shape 

■  NirC  1^2  ^  ^2) 

has  an  output  shape  of 


^out 


_ (2DT/N7rc)  [l/(t^  +  T^)] _ 

1  -  {1  -  exp[-ffr)oWL]}  exp  {-Da  [1  +  (2/7r)  tan‘^t/T)]) 


as  shown  in  Fig.  1-46. 


Fig.  1-46.  Photon  pulse  shape  at  output 
of  amplifier  normalized  to  input  for  a 
Lorentzian-shaped  pulse.  [  Reprinted  with 
permission  from  E.L.  Steele,  Optical  La¬ 
sers  in  Electronics  (Wiley,  New  York, 
1968),  p.l58.] 


A  square  pulse  going  through  the  medium  experiences  a  much  more  severe  sharpening  at 
the  front  edge; 

tU_^(t)-U_,(t-T^)] 

(2^yN)[U_^(t)-U_^(t-T^)] 

%ut^^^  “  1  -  {  1  -  exp  exp  [-2o-c4»^t] 

as  exhibited  in  Fig.  1-47. 

Independent  of  the  particular  waveform,  we  can  estimate  the  total  energy  that  a  pulse  of  a 
given  input  energy  will  extract  from  the  medium,  with  a  knowledge  of  the  small  signal  gain  and 
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Fig.  I“■47.  Ratio  of  output  photon  density 
ri|_(t)  to  input  photon  density  rj^  vs  time  in 
units  of  pulse  width  t  for  an  initially 
square  pulse  traversing  an  amplifying  me¬ 
dium  of  length  L=  10  cm,  initial  differ¬ 
ence  ofexcited  and  ground  state  densities 
Aq  =  8  X  10^^  cm “2,  and  resonance  cross 
section  a  =  2,5  X  10“2^cm2,  Total  num¬ 
ber  of  photons  per  unit  area  in  pulse  is 
q  =  4  X  10^  ^  cm"^  for  cur\'e  a,  and  q  = 
4X10^9cm"'2  for  curve  b.  [Reprinted 
with  permission  from  L.M.  Frantz  and  J.S. 
Nodvik,  J.  Appl.  Phys.^,  2348  (1963).] 


amplifier  dimensions.  Thus,  we  can  write  for  the  instantaneous  power  : 

(  /r  -W.  (t)A^  il 

Pout  +  <L/c)]  =  p.„(t)  jcy  [g^  -  (G,  -  1)  °j| 

where  the  stored  energy  in  the  medium  is 


W  =  ^  hi^AN 
o  2 


AN  is  the  initial  population  inversion,  =  e  ^  is  the  initial  small  signal  gain,  and 

w,^(t)  =  p.  (t)dt 

V/  _  00 


in 


From  the  power  expression,  we  can  estimate  the  expected  length  necessary  for  a  given  system 
and  an  effective  amplifier  cross  section.  It  can  be  shown  that  the  gain  and  energy  satisfy  the 


relation^^ 


. .  ,  W  Vw 

G.  °G.  =G. 


W.  /W 

O  O  O 


or 


W 

W  .  -  W.  .  /sinh-~-v~^\ 

L  =  out  in  ^  I _ ht^ANA  ^ 


hj^ANA  g 


g  W.  i 

,  .  ,  ^o  in  / 

h^;ANA  ! 
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We  have  used  a  program  implemented  by  S.  Marcus 
(see  Sec.  C-l-b)  to  calculate  the  inversion  for  a 
3:2:1  mixture  pumped  by  a  preionizing  e-beam  over 
a  range  of  primamy  currents,  and  a  range  of  total 
pressure  holding  E/P  constant.  Table  1-2  gives 
estimates  of  the  amplifier  dimensions  one  might 
use  for  a  multistage  amplifier  system  to  bring  the 
individual  pulses  from  2X10  ^to2J  total  energy 
(^nanosecond  pulses  and  lOmA/cm  primary  cur¬ 
rent;  the  length  restriction  is  somewhat  relaxed 
for  pulses  long  enough  to  tap  the  entire  energy 
reservoir  of  the  CO^-N^  vibrational  states).  In 
practice,  several  stages  would  be  used,  the  cross 
section  increasing  as  the  energy  in  the  pulse 
increased. 

At  this  juncture,  it  might  be  suggested  that  a 
solution  to  part  of  the  restrictions  is  to  run  the 
amplifier  system  unsaturated,  thus  eliminating  temporal  pulse  distortion  (the  unsaturated  re¬ 
quirement  being  more  easily  achieved  for  the  longer  FM  pulses).  If  the  duty  cycle  of  the  ampli¬ 
fier  is  small,  then  it  is  necessary  to  avoid  pumping  the  medium  as  hard  with  excitation  energy 
because  some  fraction  (--TS  percent)  will  always  go  into  medium  heating  and  the  remainder  will 
be  dissipated  by  this  means  also,  if  not  extracted  as  optical  energy.  Since  a  gas  (translational) 
temperature  rise  of  about  300°K  completely  shuts  off  the  gain  mechanism,  and  this  is  one  of  the 
two  significant  sources  of  pulse  distortion,  techniques  to  minimize  the  heating  over  the  range  of 
operational  modes  must  be  utilized.  For  the  long  coherence  time,  it  should  be  possible  to  re¬ 
duce  the  primary  pumping  (or  secondary  accelerating  field)  so  that  the  gain  buildup-to-optimum 
time  is  extended  to  about  the  time  that  the  next  pulse  arrives.  Thus,  medium  heating  is  reduced 
to  preserve  medium  quality  for  longer  periods.  Sufficiently  fast  flow  rates  would  aid  in  this  ef¬ 
fort.  Calculations  of  the  pumping  and  flow  requirements  are  presently  being  made.  A  serious 
source  of  medium  distortion  is  the  pressure  shock  wave  which  emanates  from  the  cathode.  One 
hypothesis  is  that  this  develops  from  sharp  temperature  differences  due  to  ion  acceleration  and 
then  ion-neutral  and  ion-cathode  collisions  in  the  cathode  fall  region.  Within  the  time  for  the 
wave  to  travel  at  the  speed  of  sound  transversely  across  the  amplifying  medium,  significant 
index-of-refraction  changes  can  occur  and  distort  the  medium.  Basic  research  is  needed  to 
study  this  effect  and  means  of  diminishing  or  eliminating  the  shock  during  the  operation  time  of 


TABLE  i-2 

SHORT-PULSE  POWER  AMPLIFIER  LENGTH 

A  (cm^) 

L  (m) 

0.01 

322 

0.03 

107 

0.10 

32 

0.30 

11 

o 

o 

3.3 

3.00 

1.2 

10.00 

0.45 

30.00  ' 

0.24 

100.00 

0.17 

the  amplifier. 

In  the  FM  pulse  sequence,  another  potential  problem  that  was  investigated  was  the  artificial 

frequency  chirping  due  to  a  strongly  varying  index  of  refraction  (anomalous  dispersion)  in  the 

44  45 

frequency  range  over  which  the  pulse  changes.  The  medium  index  of  the  CO^  is  given  by  ^ 


??(r)  =  1 


£o  iAv){v^-p) 

^  {Ap)^  [1  +  (I/Ig)]  +  ip^  -  pf- 
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where  the  first  term  is  the  VHF  limit  of  t?,  the  second  term  represents  the  significant  disper¬ 
sion  in  the  gain  region,  and  the  third  term  represents  the  contribution  from  all  other  higher  fre¬ 
quency  molecular  and  electronic  absorption  bands.  In  this  expression,  g^  is  the  line  center 
power  gain,  p  is  the  frequency  of  the  radiation,  is  the  line  center  gain  frequency,  is  the 

gain  half- bandwidth  (at  half-power),  and  c  is  the  speed  of  light.  The  second  term  could  con¬ 
ceivably  contribute  significantly  to  index  and,  hence,  wavenumber  changes  across  our  chirped 
pulse,  producing  a  phase  velocity  dependent  on  frequency,  and  thus  artificial  chirping.  Imagine 
that  we  send  in  a  pulse  whose  envelope  is  a  square  wave  of  duration  T  and  whose  fine  structure 
is  given  by 

I(t)  oc  cos  {Ztt  [p^  +  ^  (t  -  y)]  t} 

then  the  approximate  fine  structure  frequency  distribution  is  given  in  Fig.  1-48  vs  the  anomalous 
dispersion  of  the  medium  in  that  region. 

|l8-5-4023| 


Fig.  1-48.  Frequency  content  of  chirped  pulse  in  region  of  resonant  gain  and  anomalous  dispersion. 
If  we  ask  the  change  in  chirp,  we  find  that  a  chirp  of  6p/T  without  dispersion  becomes 

5p  (A  5 1  > 

-  t  t  ^ 

with  dispersion.  So  the  condition  for  the  change  in  chirp  to  be  negligible  is  for 

6j^6t  =  6p  «  1 

c 
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For  a  chirp  6p  =  2Ap  =  500MHz,  and  a  length  of  1  m, 

-15 

6pdt  =  10 
which  is  indeed  negligible. 

Power  saturation  of  the  medium  is  seen  to  simply  smooth  out  the  local  dispersion  effect. 
Saturation  will  also  cause  some  bending  of  the  wavefront  at  a  fixed  transverse  cross  section  of 
the  pulse  if  the  intensity  strongly  varies  across  that  cross  section.  It  is  clear  that  a  power  dif¬ 
ference  twice  saturation  across  the  transverse  direction  will  produce  a  phase  shift  of  half  the 
magnitude  at  the  leading  edge  {6p  =  Av)  across  the  beam.  The  index  change  due  to  dispersion, 
if  we  assume  dv  =  Av  and  g  =  0.03  cm  ^  =  13dB/m,  is  then 

Ar)(v^^Av)  =-^ 
which  corresponds  to  a  phase  shift  of 

27rAnz  =  - -  =  one  wavelength/17, 3  dB  of  gain 

(i.e.,  if  we  use  the  full  bandwidth  of  the  amplifier  for  chirping).  If  we  reduce  the  amplifier  band¬ 
width  used,  the  phase  shift  goes  down  less  than  linearly. 

H.  Granek 
R.  J.  Carbone 

2.  The  Possibility  of  Generating  Laser  Plasmas  by  Photoionization 

a.  Theoretical  Investigation 

The  possibility  of  producing  large  volume  uniform  plasmas,  via  photoionization,  with  elec- 
tron  densities  on  the  order  of  10  /cc  and  greater  in  high-pressure  gas  laser  media  is  investi¬ 
gated  as  an  alternative  to  the  e-beam  ionization  technique.  As  in  the  e-beam-excited  lasers,  the 
electrons  generated  by  photoionization  would  be  heated  in  an  externally  applied  electric  field  to 
the  optimum  temperature  for  excitation  of  laser  states. 

Two  major  problems  are  encountered.  The  first  relates  to  the  magnitude  of  the  photoioniza¬ 
tion  cross  section  a.  For  most  atoms  and  molecules,  a  is  on  the  order  of  1  Mbn  =  10~^^  cm^ 

(1  megabarn)  for  photon  energies  above  the  ionization  limit.  At  elevated  ionizing  gas  pressures, 
a  large  value  of  o*  limits  the  incident  photons'  penetration  depth  to  a  short  distance.  For  instance, 
at  1  atm  the  gas  density  N  is  approximately  3X10  /cm  ,  and  for  a  »  1  Mbn  the  penetration 
depth  L  =  1/N(t  is  «0.  3  mm.  In  this  case,  an  enormous  electron  density  is  produced  in  a  thin 
layer  of  the  gas,  which  may  be  of  use  in  special  applications  but  not  for  the  production  of  large- 
volume  plasmas. 

The  second  problem  to  be  considered  is  that  the  ionization  potentials  of  molecular  and  atomic 
species,  which  are  known  to  be  suitable  for  high-pressure  gas  lasers,  are  generally  fairly  large. 
In  these  cases,  the  corresponding  photon  wavelengths  at  energies  above  the  ionization  limit  lie 
deep  in  the  vacuum  UV,  where  transparent  windows  are  not  readily  available  and  where  absorp¬ 
tion  processes  other  than  photoionization  may  become  important  in  nonionizing  species  in  the 
laser  medium. 
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Two  methods  are  proposed  to  overcome  these  problems.  The  first  is  to  seed  the  high- 
pressure  laser  medium  with  a  low-partial-pressure  gas  having  a  low- ionization  potential  to  allow 
ionization  by  photons  at  wavelengths  larger  than  UV  cutoff  limits  of  available  windows,  and  also 
longer  than  the  strong  UV  absorption  bands  of  the  high-pressure  host  laser  medium.  The  low 
partial  pressure  allows  deep  penetration  of  the  ionizing  photon,  and  hence  uniform  production  of 
photoelectrons.  An  example  of  this  approach  is  a  1-atm  CO^-He-N^  laser  seeded  with  0.1-torr 
cesium..  The  threshold  for  photoionization  of  cesium  is  3250  A,  well  above  the  wavelengths  for 

o  O  o  4^ 

the  onset  of  absorption  in  CO,  (1700  A),  N  (1450A),  and  He  (500A).  Calculation  shows  that, 
upon  irradiations  by  a  flashlamp  emitting  1  MW/cm  between  3250  A  and  26 00  A,  electron  density 

JO  O 

in  excess  of  10  /cm  may  be  produced.  Seeding  with  cesium  may  also  be  applicable  in  the 
high-pressure  CO  laser.  By  generating  the  fourth  harmonic  of  a  Nd:  YAG  laser,  it  is  possible 
to  produce  1  J  of  UV  in  10  nsec  at  2662  A.  This  wavelength  is  near  the  peak  of  the  photoionization 
cross  section  of  rubidium.  Using  this  intense  source  to  ionize  rubidium  may  produce  electron 
densities  of  5  X  There  are  a  host  of  organic  molecules  which  have  vapor  pressure  on 

the  order  of  1  torr  at  room  temperature  and  also  have  ionization  potentials  low  enough  so  that 
UV  of  wavelengths  larger  than  1700  A  may  produce  photoionization.  Among  these  are  tripentz- 
lamine,  tripropylamine,  and  tributylamine.  The  list  of  organic  compounds  with  low  ionization 
limits  is  so  extensive  that  a  search  is  expected  to  reveal  several  that  combine  the  necessary 
characteristics  of  high  vapor  pressure,  large  cross  section,  transparency  at  the  laser  wave¬ 
length,  and  low  ionization  potential. 

The  second  method  takes  advantage  of  the  inherently  large  photoionization  cross  section  of 
most  gases  and  attempts  to  induce  photoionization  by  means  of  a  two-step  (or  multiple- step)  exci¬ 
tation  where  two  (or  several)  incident  photons  are  absorbed  by  the  high-pressure  gas  to  produce 
a  single  photoelectron.  In  a  two-step  process,  the  photon  energies  can  be  as  low  as  about  half 
the  ionization  limit,  and  hence  the  window  transmission  problems  become  less  severe.  Further¬ 
more,  in  a  multiple- step  process,  the  reduced  size  of  the  effective  photoionization  cross  section 
would  allow  deep  penetration  of  the  incident  photons  and  production  of  electrons  over  a  large 
volume.  This  method  may  be  applicable  to  the  photo  ionization  of  one  of  the  molecular  or  atomic 
species  of  the  gas  laser  medium,  or  it  may  be  applied  to  multiple- step  photo  ionization  of  an 
additional  gas  mixed  in  with  the  last  gas  laser  medium.  One  attractive  situation  is  found  in  CO 
where  the  lowest  lying  excited  electronic  states  are  triplets  lying  approximately  half-way  to  the 
ionization  limit.  The  transition  to  the  ground  state  is  only  weakly  allowed  so  the  absorption  cross 
section  is  small,  permitting  the  use  of  CO  at  1-atm  pressure.  The  intermediate  excited  states 
are  nearly  metastable  with  a  long  radiative  lifetime  so  that  large  intermediate  state  densities  may 

build  up  and  give  efficient  photoionization  in  the  second  step.  Calculation^^  shows  that,  with 

2  o  o  13/3 

1-MW/cm  illumination  between  1,500  A  and  2000  A,  electron  densities  on  the  order  of  1  0  /cm 

may  be  produced  in  1  atm  of  CO.  j  g  Levine 

A.  Javan 


b.  Experimental  Investigation 

An  ionization  cell  has  been  constructed  that  contains  a  small  linear  Xe  flashlamp  and  two 
flat  parallel  electrodes.  The  lamp  is  oriented  parallel  to  the  electrodes  in  their  midplane  and 
outside  of  the  uniform  field  region  which  exists  when  a  voltage  is  applied,  so  that  the  volume 
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between  the  electrodes  is  illuminated.  All  metal  corners  and  edges  are  removed  to  minimize 
ionization  due  to  corona  discharge.  The  cell  may  be  filled  with  various  buffer  gases  and  low 
ionization  potential  vapors  to  be  irradiated  with  an  intense  pulse  of  UV;  the  electrons  formed  by 
photoionization  are  drawn  off  by  the  electric  field  between  the  electrodes,  and  the  current  is 
measured  as  the  voltage  drop  across  a  sensing  resistor  in  series  with  the  high-voltage  supply. 
Because  the  drift  velocity  of  electrons  in  the  buffer  gas  (nitrogen)  is  known  as  a  function  of  elec¬ 
tric  field  divided  by  pressure  E/P,  the  electron  density  may  be  determined  from  the  current 
measurements  from  the  relation  j  =  nevj^  where 

3  =  current  density  =  current/area  of  electrode 
e  =  electronic  charge 

Vj^  =  drift  velocity  of  electrons  in  buffer  gas. 

The  Xe  lamp  was  supplied  by  the  Xenon  Corporation.  It  is  powered  by  a  Xenon  Corporation. 
Micropulser  and  can  be  fired  at  voltages  between  13  and  20  kV;  the  maximum  voltage  represents 
an  input  energy  of  20  J.  The  lamp  arc  length  is  8  inches  (the  electrodes  are  also  8  inches  long). 

o 

A  suprasil  envelope  was  used  to  permit  passage  of  UV  down  to  1650  A,  although  it  is  expected  that 

o 

relatively  little  energy  is  emitted  at  that  wavelength;  the  output  peaks  around  2600  A.  The  light 
pulse  has  a  full  width  at  half-height  of  1  p.sec. 

When  the  cell  is  filled  with  to  5  00  torr  and  a  potential  of  typically  lOkV  is  applied  to  the 

electrodes,  a  current  pulse  in  the  sensing  resistor  is  seen  which  closely  follows  the  light  pulse 

as  simultaneously  monitored  with  a  photomultiplier.  The  peak  current  corresponds  to  an  elec- 

tron  density  of  2  X  10  /cm  .  This  current  is  interpreted  as  the  result  of  photoemission  from  the 

-9 

aluminum  electrodes;  the  time  delay  of  less  than  10  sec  between  photon  absorption  and  electron 
emission  is  too  small  to  be  resolved. 

When  0.2-torr  tripropylamine  is  added  to  the  N^,  a  dramatic  two-order-of-magnitude  in¬ 
crease  in  the  current  flow  is  seen  upon  flashlamp  ignition.  The  current  is  found  to  vary  linearly 
with  the  electrode  voltage  between  the  limits  corresponding  to  E/P  =  0  to  E/P  =  10.  This  implies 
that  a  peak  electron  density  on  the  order  of 
10^  Vcm^  is  produced  independent  of  the  elec¬ 
tric  field.  In  turn,  this  is  strong  evidence  that 
the  enhancement  in  electron  production  is  due 
to  photoionization  and  not  to  ionization  of  tri¬ 
propylamine  by  electron  collision.  If  the  lat¬ 
ter  were  the  case,  then  an  increase  in  electron 
density  would  be  expected  as  E/P  (and  hence 
the  average  electron  energy)  is  increased; 
this,  in  fact,  is  not  observed. 

The  time  development  of  the  current  pulse 
is  markedly  different  from  the  case  of  pure 
N^.  As  seen  in  Fig.  1-49,  the  electron  density 
starts  to  grow  only  after  the  illumination  has 
reached  its  peak  and  continues  to  grow  until  the  light  pulse  terminates.  Thereafter,  the  electron 
density  decays  -  rapidly  at  first,  and  then  with  a  long  exponential  tail.  These  observations  are 
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Fig.  I“49,  Upper  trace  shows  light  pulse 
from  flashlamp  measured  by  1 P28  photo¬ 
multiplier  tube.  Lower  trace  is  current  flow 
between  electrodes  due  to  photoionization. 
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tentatively  explained  as  follows.  Ionization  takes  place  by  a  two-step  process.  This  ionization 
potential  is  about  7.5  V  so  that  two  2600-A  photons  have  sufficient  energy  to  liberate  an  electron. 
The  first  photon  populates  intermediate  states  half-way  to  the  ionization  limit;  the  second  photon 
then  raises  the  excited  electron  into  the  continuum.  The  buildup  of  intermediate  state  popula¬ 
tion  may  account  for  the  l-psec  delay  in  the  growth  of  the  electron  density.  Initial  rapid  decay 
of  the  plasma  may  be  due  to  recombination  which  proceeds  at  a  rate  proportional  to  n  .  When 
n  falls  to  a  sufficiently  low  value,  attachment  and  diffusion  become  the  dominant  loss  mechanisms 
which  proceed  at  a  rate  proportional  to  n  and  give  rise  to  an  exponential  decay. 

The  peak  electron  density  has  been  measured  as  a  function  of  flashlamp  power.  Preliminary 
results  show  that  for  relatively  low  lamp  energy,  the  electron  density  increases  approximately 
as  the  square  of  the  lamp  energy;  for  high  lamp  energy,  electron  density  increases  more  slowly  - 
approximately  linearly  with  lamp  energy.  This  behavior  is  consistent  with  two-step  excitation 
for  which  ionization  is  proportional  to  lamp  energy  squared.  At  low  lamp  energy,  and  hence 
low  electron  density,  electron  loss  varies  as  n;  thus,  the  equilibrium  electron  density  varies  as 
lamp  energy  squared.  At  high  lamp  energy,  and  hence  high  electron  density,  electron  loss  varies 
as  n^;  thus,  the  equilibrium  electron  density  varies  linearly  with  lamp  energy. 

Experimental  studies  are  being  pursued  to: 


(1)  Determine  quantitatively  the  electron  destruction  mechanisms  in 
tripropylamine, 

(2)  Study  the  photoionization  of  other  organic  molecules, 

(3)  Increase  plasma  densities  by  using  more  energetic  lamps, 

(4)  Produce  plasmas  via  photo  ionization- of  cesium  vapor, 

(5)  Produce  a  low-pressure  and  ultimately  an  atmospheric-pressure 
CO^-N^-He  photoionization  excited  laser. 

J.  S.  Levine 


3.  Other  Device  Work 


a.  The  Feasibility  of  a  10-atm  CO^-N^-He  Laser 

The  variation  in  10.6-fjLm  small-signal  gain  of  CO^-N^-He  e-beam  excited  mixtures  as  a 
function  of  pressure  is  considered.  The  well-known  formula  for  small  signal  gain  is: 


where 


a{u  =  V  ) 
o 


-  if  '*.) 


a(v  = 


^spont 


gain/centimeter  at  line  center 

population  density  in  upper  level:  g^  =  degeneracy 
population  density  in  lower  level:  g^  =  degeneracy 
spontaneous  transition  time  from  2  -►  1 
half-width  at  half-height  of  Lorentzian  gain  profile. 
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In  the  case  of  CW  pumping,  the  steady- state  population  density  N.  in  the  absence  of  a  satu¬ 
rating  optical  frequency  field  is  the  product  of  the  excitation  rate  R.  and  the  level  lifetime 
For  direct  electron  excitation  of  the  CO^  levels,  R.  is  proportional  to  the  electron  density,  the 
CO^  ground-state  density,  and  the  electron  collision  cross  section 


R. 

1 


CO2 


A  similar  relation  holds  for  the  electronic  excitation  of  nitrogen  vibrational  levels.  In  both 
cases,  the  rate  of  excitation  is  directly  proportional  to  pressure.  The  rate  of  collisional  trans¬ 
fer  of  excitation  from  nitrogen  (v  =  1)  to  CO^  (001)  is  proportional  to  the  excited-state  nitrogen 
density  =  1)  and  to  the  ground-state  CO^  density 

P  CO.  CO 


Since  N.  =  1)  is  proportional  to  N  {v  =  0),  the  collisional  transfer  rate  should  increase  with 
the  second  power  of  pressure.  At  lOatm,  the  10-|jim  transitions  are  totally  pressure-broadened; 
so,  Further,  level  lifetimes  are  inversely  proportional  to  pressure. 

So  far  it  appears  that,  in  the  CW  limit,  01  should  vary  inversely  with  pressure  with  respect 
to  direct  electron  excitation,  and  be  constant  in  pressure  with  respect  to  collisional  transfer  be¬ 
tween  N.  and  CO..  However,  in  the  case  of  e-beam  excitation,  n  is  proportional  to  the  square 
Z  Z  “ 

root  of  the  energy  loss  of  the  ionizing  beam  to  the  gas  per  centimeter  (assuming  recombination 

is  the  dominant  mechanism  for  ion  loss).  This  stopping  power  is,  in  turn,  proportional  to  the 
gas  density  or  n^  ~  n/ pressure.  Therefore,  if  sufficient  primary  electron  energy  is  available 
to  penetrate  the  medium,  the  gain  of  a  high-pressure  (collisionally  broadened)  CO^  laser  should 
exhibit  a  pressure  variation  that  is  the  sum  of  two  terms.  One  term,  due  to  direct  electron  ex¬ 
citation,  varies  as  (pressure  )~^/^;  the  other,  due  to  collisional  transfer  from  to  CO^,  should 
increase  as  (pressure 

The  parameters  appearing  in  the  gain  equation  are  evaluated  for  an  8:1:1  mixture  of  }ie:CO^:N^ 
at  10-atm  pressure  assuming  longitudinal  e-beam  excitation  by  a  10  A,  3-nsec  e-beam  source 
(Febatron  706).  Peak  gain  for  the  line  is  estimated  as  2. 2  percent/cm. 

J.  S.  Levine 


b.  Construction  of  a  10-atm  CO^-N^-He  Laser 

A  laser  has  been  constructed  based  on  the  calculations  outlined  in  the  preceding  section.  Di¬ 
mensions  of  the  active  volume  are  1X1X20  cm.  Electron-beam  excitation  is  in  the  axial  direc¬ 
tion,  parallel  to  the  optical  axis.  The  e-beam  is  focused  by  means  of  a  magnetic  lens  placed 
between  the  Febatron  706  source  and  the  2-mil  titanium  electron  window  on  the  laser.  The  mir¬ 
ror  nearest  the  e-beam  source  is  formed  by  a  thin  aluminized  mylar  film  stretched  taut  over  an 
optically  flat  pyrex  annular  ring.  This  mirror  permits  free  passage  of  the  ionizing  beam  while 
providing  96-percent  reflectivity  at  10.6  }xm.  The  far  end  of  the  laser  is  terminated  by  a  KCl 
Brewster  angle  window.  The  optical  cavity  is  completed  by  a  curved,  partially  transmitting  ger- 
tnanium  mirror. 

Further  development  of  this  system  awaits  the  availability  of  the  Febatron  source. 

J.  S.  Levine 
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c.  I^mb  Dip  in  Sealed-Off  CO  Lasers 

We  have  completed  measuring  and  evaluating  data  on  the  Lamb  dip  for  16  different  rotational- 
vibrational  transitions  in  sealed-off  CO  lasers  operating  at  a  total  pressure  of  1.5  torr.  The 
experimental  data  of  power  output  vs  frequency  were  matched  to  theoretical  curves  which  included 
relaxation  across  the  velocity  profile  as  a  major  factor  limiting  the  depth  of  the  Lamb  dip.  As  a 
result  of  these  investigations,  the  following  parameters  were  accurately  determined:  the  inhomo¬ 
geneous  and  homogeneous  linewidths,  the  relaxation  rate  out  of  the  lasing  level,  and  the  cross - 
relaxation  rate.  We  found  that  the  cross-relaxation  rates  were  of  the  same  order  in  magnitude 
as  the  dephasing  time.  A  full  report  detailing  our  work  is  nearing  completion, 

C.  Freed 
H.  A.  Haus 
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II.  OPTICAL  MEASUREMENTS  AND  INSTRUMENTATION 


A.  ZnS  UP -CONVERSION 

Additional  parameters  relevant  to  the  up-conversion  process  in  ZnS  have  been  measured> 
Whereas  previous  measurements  have  been  primarily  concerned  with  the  effects  of  DC  IR  radia¬ 
tion,  AC  effects  have  also  been  measured  now  using  chopped  IR  radiation  as  a  source  and  con¬ 
stant  UV  illumination  for  excitation.  Under  these  conditions,  an  AC  visible  output  is  observed 
on  a  DC  background.  At  frequencies  on  the  order  of  typical  TV  scanning  rates,  the  following 
results  were  obtained; 

(1 )  Peak-to-peak  AC  visible  output  decreases  about  50  percent  in  going 
from  1  to  64  Hz; 

(2)  AC  visible  output  increases  with  increasing  UV  intensity; 

(3)  AC  to  DC  background  ratio  decreases  with  increasing  UV  intensity; 

(4)  AC  signal-to-noise  ratio  (SNR)  increases  somewhat  with  increasing 
UV  intensity; 

(5)  Ac  SNR  is  1  to  2  orders  of  magnitude  worse  than  the  SNR  obtained  in 
previous  measurements  where  the  phosphor  is  exposed  to  DC  IR  radia¬ 
tion  after  (not  during)  UV  excitation.  This  is  due  to  the  large  back¬ 
ground  noise  caused  by  constant  UV  illumination  in  the  AC  case. 

These  results,  and  others,  simply  confirm  the  fact  that  in  these  materials  it  takes  a  finite  time 
(several  seconds)  to  excite  electrons  to  states  which  will  produce  the  maximum  stimulation 
effect.  Under  AC  operation,  even  at  low  frequencies,  there  is  not  enough  time  for  significant 
UV  excitation  between  IR  pulses.  In  other  measurements,  we  have  determined  that  the  spectral 
character  of  the  visible  output  is  dependent  upon  the  IR  stimulation  wavelength.  Our  results 

o 

show  that  10-|im  radiation  produces  a  band  of  visible  radiation  centered  near  5300  A,  while  2-|jLm 

o 

radiation  produces  a  band  centered  nearer  5400  A. 

Presently,  we  are  constructing  apparatus  to  determine  the  S/N  characteristics  of  the  phos¬ 
phor  screen  resolution  elements  under  uniform  IR  illumination.  We  expect  to  determine  the 
degree  of  nonuniformity  in  response  due  to  irregularities  in  particle  size,  particle  distribution 
and  layer  thickness.  The  smallest  particle  size  range  is  1  0  to  20[j.m. 

W.  J.  Scouler 


B.  IR  CAMERA  (4  inch) 

In  measurements  relevant  to  the  imaging  quality  of  the  ZnS  up-converter  phosphors,  we 
have  used  the  cooled,  4- inch-diameter,  f/l.O,  3-element-refractor  camera  shown  schematically 
in  Fig.  II- 1. 

The  camera  has  a  0,25-rad  FOV  and  is  diffraction  limited  over  the  design  wavelength  range 
of  7.5  to  12  pm.*  A  refractor  was  chosen  over  a  reflector  for  the  f/l,0  design  as  it  has  a  simpler 
and  more  compact  configuration. 


*  Optical  design  by  David  Grey  Associates,  Waltham,  Mass. 
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An  important  consideration  in  the  design  is  the  reduction  to  a  minimum  of  extraneous  IR 
radiation  at  the  image  plane.  This  was  achieved  by  cooling  to  liquid-He  temperature  and 
to  liquid- temperature  and  then  mounting  between  them  a  shutter  cooled  to  nearly  liquid-He 
temperature  and  activated  via  an  O-ring  seal  at  room  temperature.  Cooling  to  liquid- 
temperature  reduces  the  heat  load  to  the  liquid-He  section  while  reducing  emission  from 
The  cooled  shutter  prevents  any  front-end  IR  radiation  from  reaching  the  phosphor  until  the 
scene  is  to  be  viewed. 

Behind  the  phosphor  screen  are  IR  blocking  filters  which  absorb  radiation  greater  than  2  [im 
while  being  highly  transmitting  to  both  visible  and  UV  radiation.  Each  set  of  filters,  one  at 
liquid-N^  and  one  at  liquid-He  temperature,  consists  of  one  3-mm-thick  Oriel  7100  filter  and 
one  3-mm“thick  piece  of  plexiglas.  The  attenuation  of  IR  ^  2  jim  is  about  1  0  in  each  set.  The 
UV  excitation  light  and  the  visible  luminescence  from  the  phosphor  are  viewed  via  the  reflex 
camera  used  to  photograph  the  phosphor  image. 

The  4- inch- diameter  Ge  collector  lens  also  serves  as  the  vacuum  window.  A  simple 
stress  analysis  shows  that  the  atmospheric  pressure  produces  no  measurable  deterioration  of 
the  diffraction  limited  image.  is  held  in  an  adjustable  holder  to  permit  some  focusing  ability 
for  objects  closer  than  infinity.  A  bandpass  filter  directly  in  front  of  the  phosphor  can  be  used 
to  limit  the  wavelength  range  to  any  band  between  7.5  and  12  pm. 

The  reservoirs  for  the  coolants  are  1.9  and  2.7  liters,  respectively,  for  liquid  He  and  liquid 
giving  an  operating  time  in  excess  of  4|  hours.  Typical  cool  down  times  are  on  the  order  of 
3  hours. 

Mechanical  rigidity  and  thermal  isolation  of  the  three  different  temperature  sections  are 
achieved  using  specially  fabricated  ribbed  stainless  steel  tubing.  These  isolators  have  a  high 
resistance  to  twist  or  distortion  and  a  high  thermal  resistance. 

Tests  have  demonstrated  that  the  mechanical  integrity  of  the  design  is  maintained  at  cryo¬ 
genic  temperatures  and  that  its  performance  is  diffraction  limited.  This  was  determined  exper¬ 
imentally  by  replacing  the  ZnS  phosphor  screen  with  a  Ge:Cu  detector  located  behind  a  knife 
edge  in  the  focal  plane,  looking  at  a  chopped  blackbody  source  at  different  pointing  angles.  The 
variation  in  angle  from  cut- on  to  maximum  signal  is  a  measure  of  the  image  size.  From  our 
measurements  of  a  1-cm  source  at  45  m,  the  best  image  obtained  had  a  38-pm  diameter.  This 
is  to  be  compared  to  a  theoretically  calculated  value  of  34  pm,  which  is  within  experimental  limits. 

Rotation  of  the  camera  through  large  angles  produces  the  signal  vs  angle  curve  of  Fig.  II- 2. 
This  curve  was  obtained  with  a  Get  Cu  detector  located  behind  a  mask  at  the  focal  plane  defining 
an  area  on  the  detector  of  1.5  X  3  mm.  The  width  of  the  central  peak  corresponds  to  the  1.5-mm 
slit  width,  which  is  much  larger  than  the  image  size.  Thus,  the  detector  acceptance  area  for 
scattered  radiation  is  much  larger  than  is  necessary  to  view  the  direct  image.  If  the  detector 
area  were  reduced  to  just  accept  the  direct  image,  then  the  off-axis  signal  would  be  about  3  orders 
of  magnitude  lower  than  shown.  Note  that  at  angles  of  about  18°  from  the  normal,  there  are 
peaks  at  about  1/3000  of  the  signal  detected  with  the  camera  pointed  at  the  source.  Some  effort 
has  been  made  experimentally  (so  far  unsuccessfully)  to  determine  the  cause  of  these  peaks.  In 
addition,  computer  calculated  ray-traces  of  several  assumed  reflection  possibilities  have  pro¬ 
duced  no  equivalent  path  of  suspected  rays. 
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Fig.  11-2.  Detector  signal  vs  angle.  Angle  in  degrees  is  amount  blackbody  source  is  located 
off  optical  axis. 

Presently,  a  second  set  of  lenses  (L^,  and  L^)  is  being  set  up  on  an  optical  bench  at 

room  temperature  in  an  attempt  to  discover  the  path  of  the  extraneous  rays. 

D.  H.  Dickey 
T.  M.  Quist 
W.  J.  Scouler 


C.  INTERFEROMETRIC  IMAGING  TECHNIQUES 
1 .  Introduction 

The  most  serious  limitation  to  using  the  interferometric  imaging  technique  from  the  KC-135 
aircraft  in  flight  is  suspected  to  be  the  presence  of  an  image  degrading  boundary  layer.  The 
earlier  reported  retroreflector  type  measurements  showed  that  an  image  degradation  arises  be¬ 
tween  the  wing  tip  of  the  aircraft  and  the  skin.^ 

In  these  experiments,  the  measured  modulation  transfer  function  (MTF)  fell  quite  rapidly  to 
low  contrast  values,  and  there  was  some  question  raised  as  to  whether  the  MTF  actually  repre 
sented  only  the  boundary  layer  (near  the  skin  of  the  aircraft).  In  order  to  confirm  that  the  degra¬ 
dation  was  in  fact  localized  near  to  and  just  outside  the  aircraft  window,  a  new  experiment  was 
designed  and  carried  out.* 

The  boundary  layer  localized  in  the  first  25  cm  away  from  the  skin  of  the  aircraft  was  meas¬ 
ured  in  this  experiment.  In  so  doing,  it  was  shown  that  the  interferometer  performed  properly 

*  Experiment  planned  and  implemented  in  cooperation  with  D.  G.  Kocher  and  B.  W.  Bryant  in  support  of  the  LFLC  telescope  seeing 
assessment  work,  but  using  the  interferometer  as  a  sensor. 
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out  to  the  highest  spatial  frequency  (maximum  shear  value)  and  demonstrated  that  the  effects  of 
vibration,  air  currents  inside  the  aircraft  cabin  near  the  window  and  temperature  gradients 
across  the  optical  window  produced  no  significant  part  of  this  degradation. 

MTF  measurements  of  the  boundary  layer  were  made  both  at  the  front  (LFLC)  window  and 
at  the  third  window  in  the  aft  array.  The  measurements  were  then  compared  with  the  earlier 
retroreflector  results. 

2.  Experimental  Program 

A  schematic  of  the  experimental  setup  is  shown  in  Fig.  II- 3.  A  5-mW  He-Ne  CW  laser 
(wavelength  632.  8nm)  was  mounted  rigidly  on  top  of  the  interferometer.  A  high  quality  50-mm- 
aperture  laser  beam  expanding  telescope  was  used  with  the  laser  to  provide  a  collimated  and 
coherent  laser  beam.  A  mirror  and  a  beam  splitter,  mounted  to  allow  angular  adjustments 
in  each  case,  were  arranged  to  direct  the  50-mm  collimated  laser  beam  through  a  75-mm- 
diameter  optical  window  in  the  side  of  the  aircraft.  The  beam  then  propagated  through  the 


Fig.  11-3.  Schematic  of  interferometric  boundary  layer  MTF  measurement. 

air  stream  (the  boundary  layer)  to  a  gold- coated  mirror  31  mm  in  diameter  (supported  on  an 
airfoil  outside  the  aircraft)  where  it  was  reflected  back  through  the  boundary  layer  and  the  optical 
window.  The  mirror  was  located  25  cm  from  the  window  face  outside  the  aircraft.  Both  the 
window  and  mirror  were  mounted  in  an  aerodynamically  designed  airfoil  plate  unit  made  of  alumi¬ 
num,  which  can  be  fitted  into  any  of  the  aircraft  window  locations  in  place  of  the  standard  22-inch- 
diameter  optical  windows  of  the  KC-135  aircraft.*  An  aerodynamically  designed  tube  (an  airflow 


*  Airfoil  plate  design  by  R.  Stiles,  Airborne  Research  &  Engineering,  llanscom  AFB. 
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shield)  could  also  be  readily  fitted  (on  the  ground)  between  the  gold-coated  mirror  and  the  75-mm 
optical  window  to  completely  enclose  and  shield  the  light  path  outside  the  aircraft  from  the  bound¬ 
ary  layer. 

The  laser  beam  returned  through  the  optical  window,  following  the  double  pass  through  the 

boundary  layer,  and  passed  next  through  the  beam  splitter  to  a  second  smaller  24 -mm- aperture 

beam  reducing  telescope.  The  beam  was  reduced  in  diameter  to  about  an  8-mm  collimated  beam 

2 

and  then  passed  through  a  dove  prism  and  finally  to  the  corner  cube  shearing  interferometer. 

The  dove  prism  rotated  the  wavefront  around  the  optical  axis,  to  allow  measurements  for  differ¬ 
ent  diametrical  orientations  of  the  shear  axis.  The  diameter  of  the  collimated  laser  beam 
through  the  boundary  layer  was  limited  by  the  diameter  of  the  entrance  aperture  of  the  telescope 
to  a  nominal  24  mm. 

All  the  auxiliary  optical  components  arranged  in  this  experimental  setup  were  of  very  high 
optical  quality,  to  ensure  that  the  wavefront  entering  the  interferometer,  in  the  absence  of  the 
boundary  layer,  was  plane  to  within  l/lO  wave. 

A  photograph  of  the  interferometer  and  the  laser  source,  with  the  auxiliary  optical  compo¬ 
nents  in  place,  is  shown  in  Fig.  II-4(a).  The  complete  unit  was  mounted  on  a  rigid  tripod,  which 
allowed  angular  and  translational  adjustments  to  be  made  and  allowed  the  laser  beam  to  be  di¬ 
rected  through  the  75-mm-diameter  window  in  the  airfoil  plate  (a  paper  mask  in  the  photograph 
shows  the  location  of  the  airfoil  plate  and  the  two  windows).  A  photograph  of  the  airfoil  plate  in 
position  outside  the  aircraft  is  shown  in  Fig.  II-4(b).  The  laser  beam  passes  and  returns  through 
the  left  hand  window  in  Fig.  II- 4(b).  The  tripod  was  bolted  to  the  floor  to  produce  a  relatively 
rigid  mounting.  The  laser  beam  was  aimed  through  the  right  hand  window  as  seen  in  Fig.  II-4(a) 
and  a  second  window  (to  the  left  of  the  laser  beam  window)  in  the  airfoil  plate  served  simply  for 
visual  observation. 


Fig.  II -4(a).  Interferometer  at  LFLC  window 
position  in  aircraft. 


Fig.  Il■"4(b).  Airfoil  plate  from  outside 
aircraft  (airflow  shield  removed). 
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All  the  electrical  measurements  were  recorded  on  an  FM  tape  recorder  and  the  data  were 
analyzed  on  the  ground  in  the  laboratory.  An  oscilloscope  display  of  all  the  signals  provided 
continuous  real-time  monitoring  of  the  data  being  recorded  and  permitted  accurate  electro-optic 
focusing  of  the  telescope  in  flight  before  each  MTF  run. 

Both  of  the  two  airfoil  plate  windows  had  provision  for  either  a  single  or  double  (airspaced) 
optical  window  mounting  (to  circumvent  frosting  problems),  but  only  the  single  window  unit  was 
employed  in  the  experiment  reported  here. 

3.  Results 

Seven  daylight  flights  (at  about  40,000  feet  altitude)  were  carried  out.  In-flight  MTF  meas¬ 
urements  were  performed  with  the  airflow  shield  in  place  (no  boundary  layer)  and  with  it  removed 
(boundary- layer  measurement).  An  earlier  check  of  the  complete  system  was  done  in  the  labora¬ 
tory  and  in  the  aircraft  on  the  ground.  MTF  measurements  in  flight  were  made  at  the  LFLC 
window  position  (located  well  forward  of  the  wing,  7  m  from  the  aircraft  nose)  and  at  the  third 
mount  aft  array  window  position  (located  over  the  wing,  18  m  from  the  aircraft  nose).  The  inter¬ 
ferometer  was  tripod- mounted  for  the  LFLC  window  measurement,  and  was  placed  on  the  fixed 
gimballed  mount  (with  micrometer  adjustments)  for  the  third  mount  measurement. 

MTF  measurements  were  made  for  different  orientations  {2(p)  of  the  wavefront,  correspond¬ 
ing  to  the  MTF  along  a  horizontal  {Zcp  =  0),  inclined  {2cp  =  45*"  )  and  vertical  {Z(p  =  9  0® )  axis.  This 
will  suffice  to  establish  whether  there  is  any  asymmetry  in  the  wavefront  arising  from  a  non¬ 
symmetry  in  the  boundary  layer.  Each  telescope  could  be  independently  and  accurately  focused, 
and  the  focus  was  checked  before  each  MTF  run  by  maximizing  the  signal  on  the  scope  display 
at  a  shear  setting  of  approximately  s  =  1.0.  The  parameter  s  is  defined  as  a  reduced  spatial 
frequency  s  =  2XRF,  where  A  is  the  wavelength  (in  cm),  R  is  a  spatial  frequency  (in  cycles 
per  cm),  and  F  is  the  ratio  of  the  focal  length  to  aperture  diameter  of  the  optics. 

The  MTF  measurements  obtained  at  the  third  mount  aft  array  window  are  shown  in  Fig.  II- 5 
(from  3  flights).  A  similar  set  of  data  were  obtained  at  the  LFLC  position,  and  the  respective 
mean  MTF  for  each  window  position,  is  shown  compared  in  Fig.  II-6.  In  both  figures  the  diffrac¬ 
tion  limited  (clear  aperture)  MTF  curve  is  shown.  The  data  represent  a  double  pass  through 
the  boundary  layer. 

With  the  airflow  tube  in  place  (no  boundary  layer),  the  MTF  is  close  to  the  diffraction  lim¬ 
ited  curve  (curve  1,  Fig.  II-5).  When  the  tube  is  removed  (boundary  layer  present)  the  MTF  falls 
(curve  2,  Fig.  II- 5).  These  data  characterize  the  image  degradation  due  to  the  boundary  layer 
during  daylight  flights  at  40,000  feet  altitude  and  at  normal  cruising  speed.  There  is  no  signifi¬ 
cant  difference  in  data  obtained  for  different  pupil  rotations  (at  0®,  45®  and  90®)  corresponding  to 
different  shear  orientations  along  the  wavefront  diameter.  There  appears  therefore  to  be  no 
asymmetry  in  the  degradation  due  to  the  boundary  layer. 

Since  with  the  air  flow  tube  in  place  the  MTF  measured  is  virtually  diffraction  limited,  there 
cannot  be  any  significant  degradation  arising  from  extraneous  air  currents  near  the  window  inside 
the  aircraft,  nor  from  temperature  gradients  in  the  optical  window  itself. 

In  order  to  compare  these  data  with  the  theoretical  prediction  available  at  this  time,  the 
characteristics  of  the  refractive  index  (or  temperature)  fluctuations  in  the  turbulent  airstream 
must  be  evaluated.  These  in  turn  will  determine  the  random  wavefront  errors  from  which  the 
MTF  can  be  computed  for  this  case. 
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Fig.  11-5.  MTF  of  aircraft  boundary  layer  at  third 
window/aft  array. 


Fig.  11-6.  MTF  of  in-flight  aircraft  boundary  layer 
at  front  (LFLC)  and  third  window/aft  array. 


As  a  rough  estimate  of  the  wavefront  errors  in  the  present  instance  for  a  double  pass  through 
the  aircraft  boundary  layer,  a  fluctuation  for  the  optical  phase  (path  length)  of  X/7.5  is  believed 
to  be  of  the  right  order  *  This  value  will  be  used  without  detailed  discussion  to  make  a  prelim¬ 
inary  comparison  with  theory.  Barakat^  has  shown  that  the  MTF  can  be  calculated  from  random 
wavefront  errors.  Using  his  method,  the  correlation  parameter*^  was  roughly  estimated  from 
our  earlier  retroreflector  MTF  measurements.  For  (t  =  x/7.5,  the  MTF  [!/)(s)]  was  then  computed 
from  the  relationship  given  by  Barakat 

i 

\  =  exp{~k^(T^  [1  -  exp(-as^)]} 

where  ip  (s)  is  the  diffraction- limited  MTF,  s  is  the  reduced  spatial  frequency  and  k  =  2v/\. 

This  computed  curve,  for  the  conditions  spelled  out  above,  is  shown  by  the  points  on  Fig.  II-6. 
While  the  present  analysis  is  incomplete  and  should  be  regarded  as  tentative,  a  comparison  be¬ 
tween  the  measured  and  computed  curve  shows  satisfactory  agreement. 

In  order  to  compare  the  present  results  with  those  previously  obtained  in  the  wing-tip  retro- 
reflector  experiment,  in  which  data  for  effectively  a  single  pass  were  obtained  with  an  89-mm- 
diameter  beam  (with  a  central  obscuration  a  =  0.34),  a  certain  amount  of  extrapolation  is  called 
for. 

The  equivalent  single-pass  MTF  through  the  boundary  layer,  can  be  obtained  from  the  actual 
measured  double-pass  data  (Fig.  II-5)  by  assuming  that  the' relationship  between  the  MTF  {tp)  and 
the  path  length  L  through  the  boundary  layer  is  of  the  form;  L  =  k  log{^^/tp)  where  k  is  a  con¬ 
stant  and  ip^  is  the  optical  diffraction  limited  MTF.  The  equivalent  single-pass  MTF,  computed 


*  Estimated  from  sheared  interferograms  photographed  during  the  aircraft  flights. 
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in  this  way  from  the  double-pass  measured  curve,  is  next  extrapolated  from  the  24-mm-diameter 
clear  aperture  case  to  the  89 -mm- diameter  obscured  aperture  case. 

In  this  extrapolation,  it  is  assumed  that  there  is  no  inherent  aperture  size  dependence,  except 
for  scaling,  on  the  degradation.  It  is  simply  an  extrapolation  along  the  spatial  frequency  axis,  to 
correlate  this  parameter  with  the  respective  beam  diameter.  The  resulting  curve  from  this 
extrapolation  is  shown  in  Fig.  II- 7,  in  the  case  of  the  third  window/aft  array  measurement.  It  is 
compared  on  the  same  figure  with  the  MTF  measurement  made  in  the  retroreflector  experiment 
at  this  same  window  (daytime  flying).  For  the  lower  spatial  frequencies  (below  s  =  0.2)  there  is 
relative  agreement  between  the  two  curves.  The  extrapolated  curve  above  s  =  0. 2  is  not  very 
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Fig.  11-7,  MTF  of  in-flight  aircraft  boundary 
layer  at  third  window/aft  array  [comparison  of 
direct  (25-cm)  path  and  retroreflector  (19-m) 
path]  . 


0,2 


accurate,  since  it  is  determined  from  data  at  the  higher  spatial  frequencies  in  the  24-mm-beam 
experiment  where  the  SNR  is  poor.  Nevertheless,  it  provides  a  direct  comparison  of  the  in¬ 
flight  MTF  measurement  for  the  local  (25-cm)  path  and  the  retroreflector  experiment  (19  m  to 
the  wing  tip)  path  at  the  third  window/aft  array  location,  A  similar  comparison  was  made  for 
the  LtFLC  position,  with  slightly  poorer  agreement  between  the  two  curves,  due  mainly  to  the  in¬ 
fluence  of  window  astigmatism  on  the  retroreflector  MTF  curve  at  the  LFLC  position. 

4.  Conclusions 

The  boundary  layer  is  mainly  responsible  f6r  the  degraded  imaging  observed  in  flight  through 
the  two  aircraft  windows  mentioned,  the  seeing  being  slightly  better  at  the  front  window.  By  com¬ 
paring  the  local  (25-cm)  path  measurements  with  the  longer  (19-  or  27-m  retroreflector)  path 
measurements  it  is  concluded  that  the  degradation  apparently  is  confined  mostly  within  the  25-cm 
path. 

An  attempt  made  to  obtain  theoretically  computed  MTF  curves,  based  upon  present  knowledge 
(this  ultimately  has  come  from  wind  tunnel  data  to  date),  gives  reasonable  agreement  between 
theory  and  experiment  for  the  airplane.  The  assumptions  and  estimations  made  here  should  be 
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regarded  as  tentative,  until  better  substantiation  of  the  value  for  the  wavefront  errors  is  obtained. 

In  summary,  the  interferometric  system  has  been  demonstrated  to  be  capable  of  MTF  meas¬ 
urements  over  its  complete  range  (to  maximum  shear)  on  board  an  aircraft  in  flight,  and  has 
established  the  existence  of  image  degradation  due  to  a  boundary  layer. 

D.  Kelsall 


D.  AIRBORNE  IMAGE  STUDIES 

1.  Measurements  of  the  Tropopause  Effect 

In  the  fall  of  1971,  four  flights  were  made  with  the  Lincoln  Laboratory  KC-135  in  the  neighbor 
hood  of  Christchurch,  New  Zealand,  in  order  to  obtain  high  resolution  stellar  photographs  both 
above  and  below  the  atmospheric  tropopause.  The  purpose  was  to  determine  whether  any  effect 
could  be  observed  on  the  images  caused  by  the  turbulent  air  in  the  boundary  region  and  to  relate 
the  measurements,  if  possible,  to  proposed  profiles.  D.  Kelsall  has  reported  the  results  of 
MTF  measurements  which  were  made  on  the  same  flights.  He  noted  no  significant  trend  in  the 
MTF  above  and  below  the  tropopause,  nor  any  elevation  angle  dependence  (Ref.  5,  p.  47). 

A  30- cm- diameter  Schmidt- Cassegrainian  diffraction- limited  telescope  was  used  to  photo¬ 
graph  the  images  on  Kodak  plus-X  film.  Images  formed  by  this  telescope  in  flight  have  been  ex¬ 
tensively  analyzed;  the  various  sources  of  image  degradation  have  been  identified  and  measured. 

In  this  series  of  experiments  we  had  hoped  to  extract  the  "tropopause  component"  T  in  the  same 
way  that  the  other  components  were  found,  namely  by  root  differences  of  squares.  Specifically, 
two  separate  experiments  were  performed.  The  first  consisted  simply  of  photographing  a  star 
from  5  kft  below  the  tropopause  and  then  from  5  kft  above;  a  complete  focus  run  was  made  at  each 
altitude.  Determining  the  average  altitude  of  the  tropopause  was  facilitated  by  balloon  soundings 
received  every  12  hours  from  several  locations  within  a  radius  of  100  to  200  miles.  The  Christ¬ 
church  meteorological  office  was  particularly  helpful  in  obtaining  the  necessary  observations  and 
forecasting  information.  The  estimated  tropopause  position  was  checked  by  recording  the  free 
stream  temperature  vs  altitude  on  the  aircraft  during  the  experiment  and  by  noting  the  altitude  at 
which  temperature  inversion  occurred.  The  second  experiment  consisted  of  photographing  a  star 
for  10  sec  every  minute  as  the  aircraft  climbed  through  the  "disturbed"  region.  This  was  re¬ 
peated  at  a  predetermined  focus  setting  on  each  flight.  The  rate  of  climb  varied,  but  averaged 
about  1000  feet /min.;  the  total  climb  was  about  10,000  feet.  The  main  advantage  of  the  second 
experiment  over  the  first  was  that  local  weather  changes  could  be  kept  to  a  minimum  during  pas¬ 
sage  through  the  tropopause. 

After  the  photographs  were  developed,  about  5  images  were  selected  at  random  under  each 
condition,  enlarged  approximately  400X,  and  their  diameters  measured  by  estimating  the  diam¬ 
eter  containing  90  percent  of  the  exposed  grains  using  calibrated  circles.  This  method  has  been 
calibrated  by  an  actual  grain  count  to  establish  the  "eye  ball"  9 0-percent- circle.  Two  sets  of 
independent  measurements  were  made,  one  by  D,  G.  Kocher  and  the  other  by  B.  W.  Bryant.  The 
results  were  averaged  to  minimize  observer  bias.  The  results  of  the  4  flights  are  shown  in 
Fig.  II- 8.  There  appears  to  be  little  if  any  dependence  of  image  size  on  the  tropopause  position 
for  any  of  the  flights.  Other  disturbances,  possibly  local  weather,  account  for  most  of  the 
changes  in  image  size.  For  example,  the  unusually  fine  images  on  flight  3  at  about  14°  elevation 
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Fig.  1 1-8.  Image  size  vs  elevation  angle,  above  and  below  the  tropopause.  Flights  out 
of  Christchurch,  New  Zealand,  September  1971. 

are  otherwise  very  difficult  to  understand.  Of  course,  the  tropopause  is  a  very  ill- defined 
region,  in  general;  the  strength  of  the  atmospheric  turbulence  there  can  vary  greatly.  Conse¬ 
quently,  one  would  certainly  not  expect  to  find  a  constant  tropopause  component.  If  we  let 
be  image  size  below  the  tropopause  and  cp ^  be  that  above,  then  the  mean  value  of  the  component 
T  =  However,  from  Fig.  II- 8,  one  finds  that  A(p}  is  negative  (T  imag¬ 

inary)  about  40  percent  of  the  time.  Only  at  the  largest  elevation  angles  do  there  seem  to  be 
consistently  better  images  above  the  tropopause  with  T  7  X  10  ^  rad,  but  this  is  probably  il¬ 
lusory  since  7  X  lo”^  is  about  the  lower  limit  of  statistical  significance.  In  fact,  if  there  were 
a  significant  T- component  on  any  of  our  flights,  one  would  have  expected  to  see  a  greater  and 
more  consistent  difference  at  low  elevation  angles  where  the  path  through  the  disturbed  region 
is  greater.  If  one  is  located  just  at  the  base  of  the  disturbed  region,  then  the  path  length  varies 
approximately  as  cosec  of  the  elevation  angle.  The  path  at  4'*  elevation  is  about  lOX  what  it  is 
at  50°  and  therefore  T(4°  )  should  be  about  3X  T(50°  ).  All  of  our  climb-out  data  were  obtained 
at  rather  low  elevations  in  order  to  gain  as  much  sensitivity  to  the  turbulent  region  as  possible. 
Figures  II-9  through  -12  summarize  these  data.  Again,  the  fluctuations  appear  to  be  independent 
of  the  position  of  the  tropopause  within  the  sensitivity  of  the  measurements.  Only  in  flight  4 
(Fig.  11-12)  where  the  elevation  angle  was  4.1°  is  there  even  a  suggestion  of  smaller  images  in 
the  stratosphere.  These  data  lead  to  T  ^  13  x  10”^-rad  as  an  upper  limit  for  the  contribution  of 
the  stratospheric  boundary- layer  turbulence  to  the  90-percent  image  size. 

Admittedly,  these  data  represent  a  small  sample  at  one  location  and  time  of  year.  However, 
similar  data  were  collected  near  Anchorage,  Alaska,  in  August  1968.  Though  not  as  extensive, 
the  results  were  very  similar  to  those  shown  in  Fig,  II- 8.  The  main  limitation  to  the  sensitivity 
of  this  method  is  the  existence  of  a  large  and  variable  "airborne  seeing"  component  (now  identified 
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Fig.  11-9.  Image  size  vs  altitude,  flight  1, 
4  September  1971 


Fig.  11-10.  Image  size  vs  altitude,  flight  2, 
6  September  1971 . 


Fig.  11-11.  Image  size  vs  altitude,  flight  3, 
9  September  1 971 . 


Fig.  11-12.  Image  size  vs  altitude,  flight  4, 
12  September  1 971 . 
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principally  with  the  aircraft  boundary  layer),  which  causes  the  images  to  be  much  larger 
(~20  X  lo”^  rad)  than  the  Airy  disk  (~4  X  lo’^  rad).  The  causes  for  its  fluctuation  are  not  known, 
but  may  be  related  to  the  air  mass  changes  in  the  upper  atmosphere.  During  relatively  quiet 
periods,  the  sensitivity  limit  is  determined  by  the  statistical  fluctuations  in  the  measurements 
(ct  1  X  lO^^rad)  of  these  "enlarged"  images.  Therefore,  the  limiting  sensitivity,  S^j^mit  ^ 

J (21  )^  -  (20)^  6.3  X  lo"*^  rad;  i.e.,  any  component  less  than  this  is  lost  in  the  "noise." 

Because  additional  noise  sources  are  always  present,  the  practical  limit  is  probably  closer  to 
8  or  9  X  lo”^  rad.  These  results  are  currently  being  compared  with  various  theoretical  models 
of  the  upper  atmosphere. 

2 .  Boundary  Layer 

Another  series  of  aircraft  experiments  was  planned  and  carried  out  in  March  and  April  of 
1972  which  showed  convincingly  that  the  turbulent  boundary  layer,  a  region  about  10  to  25  cm 
thick,  just  beyond  the  outer  surface  of  the  KC-135  fuselage  is  a  major  cause  of  airborne  image 
degradation.  The  measurements  were  made  using  an  interferometer,  which  was  originally  de¬ 
signed  by  D.  Kelsall  for  another  purpose,  and  which  fortunately  became  available  in  time  for  its 
inclusion  in  the  planning  of  a  boundary- layer  experiment.  The  description  of  the  experimental 
setup  and  the  salient  results,  some  of  which  are  used  here,  are  discussed  in  Sec.  II- C  of  this 
report.  Here,  we  compare  one  of  these  results  with  the  residual  B-term  (Ref.  1.,  p.  42)  in  the 
experimental  formula  for  the  90-percent  image  diameter;  this  term  was  tentatively  ascribed  to 
boundary- layer  turbulence.  Experimentally,  the  magnitude  of  B  is  approximately  18X10  rad. 
The  point  spread  function  of  the  boundary  layer  can  be  found  by  extracting  the  Fourier  transform 
of  the  single-pass  boundary- layer  MTF  measured  at  the  LFLC  window  (see  Sec.  II-C).  Figure  11-13 
illustrates  the  MTF  as  far  as  the  diffraction  limit  of  the  1-inch  probe  beam  will  allow  its  compu¬ 
tation.  The  MTF  was  computed  by  dividing  the  diffraction- limited  MTF  into  the  "measured"  MTF 


Fig.  11-13.  Single-pass  boundary  layer  MTF  with  diffraction  limit  divided  out. 
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(slightly  corrected  for  noise  at  the  high  spatial  frequencies),  which  included  both  boundary  layer 
and  diffraction.  Unfortunately,  the  MTF  rapidly  becomes  indeterminate  near  the  diffraction 
limit,  as  indicated  by  the  error  bars  which  are  rough  estimates  of  the  uncertainties.  The  curve 
was  extended  to  high-frequency  cutoff  by  an  arbitrary  linear  extrapolation,  which  is  probably  not 
unreasonable,  but  obviously  not  very  accurate. 

The  transform  was  carried  out  numerically  using  the  Cooly-Tukey  fast  Fourier  transform  al¬ 
gorithm.  The  resultant  point  spread  function  (Fourier  amplitude  squared)  is  shown  in  Fig.  11-14 
along  with  the  90  percent  width,  which  indicates  the  points  on  a  2-dimensional  curve  which  corre¬ 
spond  to  the  circle  on  a  3- dimensional  curve  of  rotation  which  enclosed  90  percent  of  the  volume. 


Fig.  11-14.  Boundary  layer  point  spread  function  with  indication  of  90-percent  image  width. 

The  abscissa  can  be  converted  to  an  angular  scale  by  dividing  by  the  focal  length  (85  mm)  of 
the  entrance  optics  to  the  interferometer.  The  90-percent  angular  width  was  computed  to  be 
17.5  X  lO'^rad,  which  is  close  enough  to  18  x  lo"^  to  leave  little  doubt  that  the  B-term  repre¬ 
sents  the  boundary- layer  component. 

A  number  of  persons  participated  in  various  aspects  of  these  experiments;  among  them  are 

D.  G.  Kocher,  D.  Kelsall,  E.  A.  Davidson  and  Lt.  R.  Stiles. 

B.  W.  Bryant 

E.  TUNABLE  INFRARED  LASER  SPECTROSCOPY  OF  SULFUR  DIOXIDE 
1 .  Introduction 

Sulfur  dioxide  (SO^)  is  one  of  the  most  common  and  harmful  of  the  atmospheric  pollutant 
gases.  With  the  advent  of  laser  techniques  for  remote  or  in  situ  monitoring  of  SO^  and  other 
gases  by  their  IR  "signatures,"^  knowledge  of  their  fundamental  spectra  becomes  imperative. 

As  we  described  earlier,"^  tunable  semiconductor  diode  lasers  are  being  developed  for  air 
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pollution  monitoring  applications  as  well  as  for  ultra -high- re  solution  IR  spectroscopy.  Doppler- 
limited  spectra  have  already  been  obtained  for  several  gases  in  different  spectral  regions:  sulfur 
hexafluoride  (SF^),  ammonia  (NH^),  and  ethylene  (C^H^)  around  10.6  jam  using  Pb^  ggSn^ 
tunable  diode  lasers;^  and  carbon  monoxide  (CO),  nitric  oxide  (NO),  and  water  vapor  in  the 

9 

4- to  6-fxm  region  using  PbS.  „Se  diode  lasers. 

10^ ^ 

In  an  earlier  report  we  showed  preliminary  diode-laser  studies  of  the  band  of  SO^, 
centered  around  8.7  pm.  The  purpose  of  the  present  report  is  to  show  agreement  between  the 
experimental  data  and  theoretical  predictions  based  on  accurate  microwave  absorption  data;  and 
thereby  to  establish  a  new,  more  accurate  value  for  the  band  center.  We  also  show  how  these 
results  can  be  used  to  predict  the  IR  spectra  of  this  gas  under  varying  conditions  of  temperature 
and  pressure,  so  that  optimum  wavelengths  for  monitoring  applications  can  be  selected. 

2.  Diode  Laser  Fabrication  and  Experimental  Procedure 

-1 

In  order  to  produce  diode  lasers  which  emit  between  1100  and  1200  cm  ,  single  crystals  of 

11 

Pbi-xSn^Te  were  grown  from  the  vapor  phase  by  a  closed-tube  process,  at  three  different 
compositions:  x  =  0.0575,  0.0700  and  0.0745.  The  crystals  were  cleaved  into  rectangular  paral¬ 
lelepipeds,  each  containing  a  p-n  junction,  of  average  size  0.084  X  0.025  X  0.020  cm  for  the  six 
diode  lasers  used  in  this  study.  Low- resistance  ohmic  contacts  were  formed,  and  the  lasers 
were  mounted  onto  the  cold-finger  of  a  liquid-He  Dewar.  Fine-tuning  of  the  wavelength  within  a 
single  longitudinal  mode  of  the  laser  was  achieved  by  changing  the  diode  current,  which  alters 
the  refractive  index  of  the  semiconductor  through  heating  of  the  junction  region. 

Each  of  the  six  lasers  exhibited  a  maximum  total  CW  power  from  the  front  face  in  the 
6  to  44-|j,W  range.  Absorption  spectroscopy  was  performed  by  transmitting  the  tunable  laser 
radiation  through  a  30- cm- long  cell  containing  the  SO^,  through  a  grating  spectrometer  and  onto 
a  Ge:  Cu  detector.  (The  purpose  of  the  spectrometer  is  to  provide  approximate  calibration  of 
the  laser  wavelength  and  to  limit  the  IR  bandpass  to  a  single  longitudinal  mode. )  The  spectra 
were  displayed  on  an  X- Y  recorder,  whose  X-axis  was  driven  by  a  signal  proportional  to  the 
diode  current  (usually  quite  linear  with  respect  to  the  wavenumber  of  the  emitted  radiation).  The 
Y-axis  was  driven  by  a  signal  proportional  to  the  transmitted  IR  power,  as  monitored  on  the  de¬ 
tector  by  a  lock-in  amplifier  synchronized  to  a  mechanical  chopper  at  the  spectrometer  entrance 
slit. 

The  SO^  was  obtained  from  the  Matheson  Company,  Inc.  of  East  Rutherford,  New  Jersey, 
with  a  quoted  purity  of  99.98%.  Gas  pressure  was  measured  with  a  capacitance  manometer  at  and 
below  1  torr,  and  with  a  Wallace -Tiernan  gauge  at  higher  pressures.  The  usual  test  sequence 
consisted  of  scans  of  SO^  at  o-,  0,1-,  1-  and  10-torr  pressure,  followed  by  NH^  for  calibration 

purposes.  Use  was  made  of  recent  precise  measurements  of  the  NH,-line  positions  in  this  re- 

12  ^ 

gion  made  available  to  us  by  Rao,  Curtis  and  Yin. 

3 .  Theory 

Because  of  the  thousands  of  relatively  strong  rotational  transitions  in  identification  of 

individual  lines  in  the  p.  vibration- rotation  band  is  somewhat  difficult.  Fortunately,  there  are 

^  13 

accurate  microwave  absorption  data  for  both  the  ground  and  the  first  excited  vibrational  states. 

For  both  vibrational  states,  the  microwave  measurements  can  be  used  to  evaluate  the  constants 
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occurring  in  a  symmetrized  expansion  of  the  rotational  Hamiltonian  in  powers  of  the  angular 
momenta  For  values  of  J  up  to  60,  the  energy  levels  as  well  as  the  transition  frequencies 

and  strengths  for  the  ^1.-3  transitions  have  been  calcu¬ 

lated.  (P,  Q  and  R  have  their  usual  meaning,  with  the  first  and  second  subscripts  referring  to 
AK^  and  AK^^,  respectively. )  In  order  to  calculate  band  strengths  from  measurements  on  a 
single  spectral  line  at  room  temperature,  the  rotational  partition  function  was  evaluated  and 
found  to  be  =  5805. 

4.  Results 

.  -1 

Figure  11-15  shows  a  1.  4-cm’^-wide  scan,  between  1156.3  and  1157.7cm  ,  for  SO2  at  pres¬ 

sures  of  0.1,  1  and  10  torr,  together  with  two  NH^  lines  for  calibration.  The  10-torr  curve 
shows  definite  pressure-broadening  and  overlap  of  closely- spaced  lines  separated  at  the  lower 
pressures.  Absorption  at  0.1  torr  is  difficult  to  observe  because  of  the  noise  level.  However, 
it  can  be  seen  for  several  of  the  strong  transitions.  Some  spectra  in  other  wavelength  regions 
where  the  laser  power  was  greater  were  considerably  less  noisy,  and  stronger  transitions  than 
those  shown  in  Fig.  11-15  could  be  easily  seen  at  pressures  below  0.1  torr.  The  individual  SO^ 
transitions  were  identified  by  comparing  their  relative  positions  and  intensities  with  theory,  and 
some  of  the  stronger  lines  in  Fig.  11-15  are  tabulated  in  Table  II- 1.  Less  than  5  percent  of  the 
lines  strong  enough  to  be  measured  by  the  20  separate  spectral  regions  covered  with  the  tunable 
diode  lasers  could  not  be  explained  by  theory.  In  order  of  decreasing  probability,  we  estimate 
these  lines  to  be  caused  by:  (a)  "hot"  bands,  the  state  being  8-percent  occupied  at  room  tem¬ 
perature,  (b)  other  isotopes  and  (c)  weaker  transitions  not  included  in  the  calculation. 


DIODE  CURRENT  (A) 


Fig.  11-15.  Diode  laser  scans  of  SO2  at  pressures  of  0,  0.1,  1  and  10  torr  at  room  temperature 
with  a  30-cm  cell .  Scan  of  NH3  at  1  torr  was  also  made  for  calibration.  Diode  laser  excita¬ 
tion  (and  tuning)  current  is  shown  by  upper  horizontal  scale.  Wavenumber  scale  represents 
theoretical  positions  of  observed  lines  (see  Table  ll-l)  using  1151,71  cm  ^  as  band  center. 
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TABLE  ll-l 

IDENTIFICATION  OF  LINES*  IN  FIGURE  11-15 


Trans  itiont 

Predicted  Wavenumber 
(cm“l ) 

Predicted  Relative 
Absorption  Coefficient 

^1,5  ^  ^2,4 

1156.335 

3.55 

^’5,17  ’  ^°6,14 

.336 

1.11 

■*  ’®2,16 

.523 

7.82 

CN 

t 

CO 

.580 

2.18 

"  ’®2,16 

.615 

2.14 

^^4,20  ""  ^^3,2] 

.780 

1.56 

’Vl4  -  1^1,13 

.841 

4.67 

’1,9  ^Vio 

.911 

5.58 

^1,3  “  ^2,2 

.979 

1.40 

^°5,15  ■"  ’\l4 

1157.056 

1.10 

^0,6  ^1,7 

.141 

4.15 

’\l2  -  ’"5,9 

.146 

1.10 

®1,5  ^  ^2,4 

.183 

2.33 

20  20 
^n,19  2,18 

.222 

6.94 

22  -►  22 
^^2,20  3,19 

.442 

7.63 

24  -►  24 

"^7,22  3,21 

.455 

7.21 

^^2,12  ■"  ’"^1,13 

.478 

3.34 

^1,7  ■*  ^2,6 

.481 

3.01 

*  Only  those  lines  for  which  the  relative  absorption  coefficients  are  greater 
than  1.0  are  included. 

t  Lower  energy  state  is  listed  first. 
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WAVENUMBER  (cm‘S 


Fig.  Il”16.  Comparison  between  measured  and  predicted  strengths  and  positions  for  some  lines  in  band 
of  SO2.  Scans  were  computer-generated  using  Lorentzian  line  profiles.  Lines  designated  as  a,  b  and  c 
represent  transitions  287  21  -29^,24'  ^8  26^^7,29  273  25  -  264^22.  respectively,  which  ap¬ 

pear  to  be  shifted  slightly  to  lower  wavenumfeer  (see  text).  Line  d  is  not  predicted  by  theory. 


PRESSURE  (torr) 


Fig.  II-IZ.  Self-broadened  and  air-broadened  dependence  of  full  width  at  half  maximum 
intensity  for  two  lines  at  room  temperature:  13 '^^^4  ]2  1163.297  cm  and 

8,  ,-8_  at  1148.810cm”’. 
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The  present  results  yield  a  value  of  1151,71  =fc  0.01  cm  for  the  i^^-band  center  frequency 
of  which  can  be  compared  with  1151.38  cm"^  obtained  by  Shelton,  Nielsen  and  Fletcher.^^ 

The  uncertainty  of  ±0.01  cm"  in  our  determination  is  due  in  part  to  possible  errors  in  the  NH^ 
spectra,  the  theoretical  calculations,  and  the  inter  comparison  between  the  NH^  and  SO^  traces. 

Figure  II-16  represents  a  detailed  comparison  between  experimental  data  from  a  diode  laser 
scan  in  the  1147”Cm"^  region  and  the  theoretical  prediction  for  lines  in  this  region.  The  curves 
marked  ’’experimental”  are  Lorentzian  fits  to  the  line  intensities,  incorporating  a  linewidth  which 
best  approximates  the  shape  of  the  experimental  curve  for  the  1-torr  pressure  at  which  the  data 
were  obtained.  (Actually,  at  l-torr  pressure  the  SO^  line  profile  is  a  Voigt  convolution  of 
Lorentzian  and  Gaussian  shapes. )  There  are  minor,  but  significant  differences  between  the  two 
curves,  as  indicated  at  a,  b  and  c,  where  transitions  of  high  J  value  appear  to  be  shifted  to 
the  left  in  the  experimental  curve  by  a  few  MHz/unit  change  in  J.  At  point  d  is  a  line  which  is 
not  predicted  theoretically  and  may  be  caused  by  one  of  the  reasons  mentioned  above. 

Figure  11-17  illustrates  the  pressure-broadening  rates  for  two  SO^  lines.  Self- broadening 
of  the  ISj  ^3  -  15_J  ^2  1163.297  cm"^  occurs  at  a  rate  of  40  ±  5  MHz/torr  for  the  full 

width  at  half- maximum  at  pressures  above  a  few  torr.  Similarly,  the  effect  of  increasing  the 
partial  pressure  of  air  on  an  SO^  line  at  fixed  pressure  is  shown  for  the  8^  .j,  g  transition  at 

1148. 81 0  cm”^,  where  the  broadening  rate  is  only  12  ±  2  MHz/torr.  Although  Fig.  11-17  does  not 
represent  linewidth  measurements  on  the  same  transition,  measurements  of  the  widths  of  sev¬ 
eral  lines  do  not  indicate  any  noticeable  variation  from  the  pressure-broadening  rates  shown 
here.  Consequently,  the  ratio  of  3.3  ±  0.7  between  self-broadening  and  air-broadening  rates  is 
significant  and  can  be  compared  to  the  5  ±  1  obtained  by  Burch,  et  ^. ,  when  dry  is  used  in¬ 
stead  of  air.  (The  presence  of  water  vapor  in  the  laboratory  atmosphere  used  for  our  air¬ 
broadening  studies  may,  in  fact,  be  the  reason  for  the  smaller  ratio. ) 


By  using  measured  absorption  constants  and  pressure-broadened  linewidths  at  10-torr  pres¬ 
sure  for  several  SO^  lines,  we  have  calculated  their  fundamental  line  strengths.  On  the  basis 

of  these  determinations,  and  using  the  theoretical  ratios  of  total  partition  function  to  line  strength, 

-20  -1  "1  2 

the  total  strength  of  the  band  is  found  to  be  358  ±20X10  cm  molecule  cm  .  This  can 
be  compared  with  the  value  of  371  ±  20  X  lo’^^cm"^  molecule"^  cm^  measured  by  Burch,  ^  ^., 
and  is  in  the  range  found  by  earlier  workers^  using  conventional  techniques.  An  effective 
dipole  moment  of  0.  086  ±  0.  003  Debye  is  found  by  a  similar  calculation. 


5.  Application  of  Theoretical  Spectra  to  Pollutant  Gas  Monitoring 

Knowledge  of  the  strengths,  wavelengths  and  widths  of  individual  vibration- rotation  lines  of 
gaseous  pollutants  can  be  used  to  predict  the  IR  spectra  corresponding  to  a  variety  of  pressure 
and  temperature  conditions.  For  the  case  in  which  SO^  in  smokestack  effluent  is  to  be  monitored, 
the  gas  is  at  atmospheric  pressure  and  can  have  wide  temperature  fluctuations.  Using  the  theo¬ 
retical  data  for  the  thousands  of  lines  in  the  8.7-|jim  band  of  SO^,  it  was  possible  with  the  aid  of 
a  digital  computer  to  display  the  absorption  spectra  for  several  temperatures.  A  segment  of 
this  band,  from  1110  to  1142  cm”^,  is  illustrated  in  Fig.  11-18  for  selected  temperatures  from 
0  to  400*0.  Below  1126  cm"^,  absorption  is  highest  for  higher  gas  temperatures,  whereas  the 
converse  is  true  above  1126  cm"^.  Consequently,  IR  absorption  of  SO^  has  the  smallest  temper¬ 
ature  dependence  around  1126  cm”^,  which  makes  this  region  particularly  appropriate  for  in- stack 
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Fig.  11-18.  Computed  spectra  for  a  portion  band  of  SO2  in  air  at  atmospheric  pressure, 

at  temperatures  from  0®  to  400®C.  A  constant  width  of  0.3  cm"^  was  used  for  all  lines. 

-1 

monitoring  of  SO^.  Moreover,  by  alternately  tuning  the  laser  from  1125.8  to  1126.8cm  and 
electronically  measuring  the  logarithm  of  the  ratio  of  the  transmitted  signals  at  the  two  wave¬ 
lengths,  the  resulting  voltage  is  directly  proportional  to  pollutant  concentration  and  is  not  affected 
by  absorption  or  scattering  processes  which  are  invariant  over  the  small  wavelength  interval. 
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F  EXTENSION  OF  TUNABLE  LASER  TECHNIQUES 
TO  AIR  POLLUTION  DETECTION 

We  have  continued  to  apply  techniques  using  tunable  infrared  lasers  to  the  problem  of  de¬ 
tection  and  measurement  of  air  pollution.  As  part  of  this  effort  we  have  considered  the  use  of 
acousto -optic  absorption  cells  in  conjunction  with  discretely  tunable  gas  lasers  and  tunable  diode 
lasers.  We  have  used  tunable  diode  lasers  for  in-situ  measurements  of  auto  exhaust,  as  well  as 
for  long  path  studies  of  background  concentrations.  The  work  described  concerns  the  gas 
ethylene  (C^H^),  although  we  have  also  studied  SF^,  NO  and  CO. 


*  Author  not  at  Lincoln  Laboratory. 
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1.  Spectrophone  Investigation 

As  an  alternate  method  to  tunable  semiconductor  laser  transmission  measurements  for 
pollutant  detection,  we  have  considered  a  technique  whereby  an  acoustic  signal,  driven  by  the 
absorption  of  IR  radiation,,  is  used  to  measure  trace  gas  concentration.  The  technique,  first 
investigated  in  the  19th  Century,^®  has  been  used  recently  to  measure  decay  times  in  gases 
and  to  detect  traces  of  methane  and  nitric  oxide  in  air  samples. 

The  acousto-optic  absorption  cell,  or  spectrophone,  can  assume  a  variety  of  forms.  The 
constructions  we  have  used  consist  of  an  inner  and  an  outer  chamber  separated  by  a  capacitance 
microphone.  The  air  sample  is  admitted  to  the  outer  chamber  and  diffuses  into  the  inner  cham¬ 
ber  by  means  of  a  small  hole.  Laser  radiation  passes  in  and  out  of  windows  attached  to  the 
inner  chamber.  If  the  laser  frequency  coincides  with  an  absorption  line  of  the  gas,  some  of  its 
energy  is  lost  to  the  gas  and  appears  as  translational  energy  causing  a  pressure  differential 
between  the  inner  and  outer  chambers.  By  chopping  the  laser  beam  and  synchronously  detecting 
the  capacitance  change  produced  by  this  pressure  differential,  we  have  measured  ethylene  con¬ 
centrations  as  low  as  28ppb.  This  measurement  used  a  0.5W  stable  CO2  laser  timed  to  the  P(14) 
line,  which  nearly  coincides  with  an  absorption  line  in  ethylene. 

Unfortunately,  at  this  low  concentration  a  signal  produced  by  absorption  in  the  windows  is 
dominant  and  must  be  subtracted  from  the  measured  signal  to  calculate  the  effect  of  the  trace 
gas.  In  this  particular  example,  the  (BaF2)  window  signal  was  2.5X  the  intensity  of  the  signal 
due  to  28ppb  of  C2H^.  Consequently,  a  high  degree  of  stability  in  laser  amplitude  and  in  con¬ 
figuration  alignment  is  required  to  complete  the  measurement. 

This  is  not  the  case  with  lower  power  lasers  where  the  sensitivity  of  the  spectrophone  deter¬ 
mines  the  lower  limit  of  detection.  Several  of  the  spectrophone  cells  we  constructed  were 
capable  of  detecting  less  than  lO’^W  of  absorbed  power  (in  the  absence  of  a  large  window  absorp¬ 
tion  signal).  With  a  5-pW  PbSnTe  diode  laser  we  were  thus  able  to  detect  an  a'cLofZXlO  , 

where  the  length  L  used  was  10  cm,  c  is  the  gas  concentration  (ppm)  and  «'  is  its  absorption 
coefficient  in  (cm- ppm)  Since  a'  is  typically  10  ^/cm  •  ppm  for  C2H^.  this  corresponds  to 
a  minimum  detectable  concentration  near  200 ppm. 

In  Fig.  11-19  we  show  a  comparison  between  the  spectrophone  output  and  the  signal  received 
by  a  Ge:Cu  detector  after  the  B-pW  PbSnTe  laser  radiation  has- passed  through  the  spectrophone. 
The  diode  radiation  was  tuned  over  a  region  near  10.6  pm,  and  the  Doppler- limited  spectrum  is 
evident  in  both  the  spectrophone  and  the  transmission  signals.  While  the  minimum  detectable 
concentration  using  the  spectrophone  was  more  than  10  ppm,  a  detection  limit  of  1  ppm  was 
possible  with  the  liquid- He -cooled  Ge;  Cu  detector  under  these  conditions.  This  reflects  the  fact 
that  the  spectrophone  NEP  was  near  10"®W,  while  the  GeCu  detector  had  a  NEP  of  approximately 
3  X  lO'^^  W. 

2.  In-Situ  Auto  Exhaust  Measurements 

The  diode  laser,  being  compact,  stable  and  easily-operated,  is  highly  suitable  for  field 
measurements.  Utilizing  this  capability  we  have  performed  fast  response  in-situ  measurements 
of  the  ethylene  concentration  in  auto  exhaust.  The  undiluted  exhaust  was  drawn  through  a  win¬ 
dowless  1.1 -m  sample  cell  and  exhausted  with  a  standard  garage  exhaust  fan.  Radiation  from  a 
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PbSnTe  diode  laser  operating  at  10.57|Jim  was  collimated,  passed  through  the  sample  cell,  and 
collected  by  a  Ge:  Cu  detector.  The  detector  output  was  amplified  with  a  logarithmic  amplifier 
and  fed  to  a  strip  recorder  whose  deflection  was  proportional  to  concentration. 

Calibration  was  accomplished  by  substituting  a  cell  with  known  ethylene  concentration.  Cali¬ 
bration  and  lack  of  interference  from  other  gases  was  further  checked  by  low  pressure  (Doppler 
regime)  studies  of  exhaust  samples  and  samples  of  known  concentration. 

A  typical  measurement  sequence  is  shown  in  Fig.  11-20.  In  Fig.  II-20(a)  the  emission 

from  a  1972  station  wagon  is  measured  after  a  "normal"  start.  The  emission  shows  two  con¬ 
centration  peaks  after  which  the  level  approaches  350  ppm.  Figure  II-20(b)  shows  a  "rich" 
start  for  which  the  accelerator  pedal  was  depressed  once  before  starting.  The  emission  exceeds 
2000 ppm  briefly  then  approaches  the  35  0  ppm  steady- state  level.  Removal  of  the  air  cleaner 
on  this  vehicle  leaned  the  fuel-air  mixture  and  reduced  the  steady  state  ethylene  concentration 
to  I60ppm.  Other  vehicles  tested  showed  somewhat  higher  concentrations.  Figure  11-21  shows 
the  ethylene  content  of  a  poorly  tuned  1964  sedan.  The  average  concentration  is  800  ppm.  The 
wide  trace  is  not  due  to  system  noise,  but  rather  to  rapid  fluctuations  in  this  car's  ethylene 
output. 

These  tests  were  designed  to  demonstrate  the  suitability  of  diode  laser  transmission  studies 
to  real-time  measurement  of  pollution  sources.  They  will  be  extended  to  include  other  gases, 
such  as  nitric  oxide  and  carbon  monoxide. 


1  ^  I - r 
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Fig.  11-21.  On-line  ^2^^  emission  test  for  untuned  1964  sedan  (system  using  pulsed 

Pb*  diode  laser  at  10:57  urn), 

U.  oo  U.  I  2 
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3.  Long-Path  Laser  Transmission  Studies 

To  increase  the  sensitivity  of  the  diode  laser  transmission  technique  to  measure  pollutant 
levels  in  ambient  air  (typically  10  to  lOOppb)  it  is  necessary  to  use  longer  path  lengths  than  are 
convenient  for  laboratory  application.  Also,  the  integrated  concentration  obtained  from  long 
path  measurements  is  often  more  suitable  than  point  by  point  results  for  regional  air  quality 
studies. 

We  have  operated  the  FbSnTe  diode  laser,  described  above,  over  a  0.5-km  path  from  Lincoln 
Laboratory  to  a  remote  reflector  and  back.  The  experimental  path,  illustrated  in  Fig.  II” 22, 
passed  over  the  Laboratory  parking  lot  and  allowed  us  to  measure  changes  in  the  ethylene  content 
during  peak  traffic  hours.  Figure  11-23  shows  an  increase  in  the  measured  concentration 

of  over  1  ppm  during  the  late  afternoon  departure  of  Laboratory  personnel.  Peaks  due  to  individ¬ 
ual  cars  (such  as  the  poorly  tuned  1964  sedan  mentioned  above)  were  also  noted. 

In  this  experiment,  the  laser  radiation  was  collimated  to  a  2  X  4- cm  beam  before  transmis¬ 
sion.  It  was  collected  by  a  larger  mirror,  used  off-axis,  and  concentrated  on  the  liquid- He- 
cooled  detector.  With  this  very  simple  arrangement,  we  were  able  to  collect  more  than  30  per¬ 
cent  of  the  transmitted  energy.  The  SNR  obtained  was  limited  by  atmospheric  fluctuations, 
rather  than  diode  power  or  detector  considerations.  As  a  result  of  these  measurements,  we 
expect  to  be  able  to  extend  the  path  length  to  many  kilometers,  without  difficulty.  By  using  the 
two-pulse  frequency  scanning  technique,  developed  for  the  EPA  smokestack  measurement  appara¬ 
tus,  we  expect  to  obtain  detection  sensitivities  near  1  ppb. 

H.  A.  Pike 
E.  D.  Hinkley 
A.  R.  Calawa 
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HL  LASER  RADAR  AND  TRACKING 


A.  INSTRUMENT  FOR  ATMOSPHERIC  MODULATION  TRANSFER 
FUNCTION  MEASUREMENTS 

1.  Introduction 

A  modulation  transfer  function  (MTF)  instrument  which  is  of  simple  design  and  easily 
transported  has  been  constructed  for  visible  atmospheric  propagation  measurements.  A  picture 
of  this  instrument  and  a  schematic  drawing  of  its  operation  are  shown  in  Figs.  III-l(a,  b). 

An  8.6-cm-diameter  Questar  telescope  is  used  to  form  an  image  of  a  distant  point  source, 
e.g.,  a  He-Ne  laser.  This  point  source  is  re-imaged  onto  the  rotating  MTF  reticle,  which  is  a 
Siemann  star  pattern  of  alternately  reflecting  and  nonreflecting  2°  segments  or  bars.  The  bar 
spacing,  that  is,  the  spatial  frequency,  is  a  function  of  radius  on  the  reticle  and  is  selected  by 
tilting  the  scan  mirror.  The  modulation  of  the  reflected  light  at  each  radius  setting  produced  by 
the  spinning  reticle  and  detected  at  the  photomultiplier  gives  a  direct  measure  of  the  MTF  at 
that  spatial  frequency. 


Fig.  MI-1.  MTF  instrument:  (a)  photograph;  (b)  schematic  diagram . 
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Fig.  1 11-2.  Instrumental  modulation 
transfer  function. 


Figure  III~2  shows  the  theoretically  calculated  MTF  for  the  instrument  along  with  a  corre¬ 
sponding  experimental  measurement.  The  particular  shape  is  due  to  the  central  obstruction  in 
the  Questar  aperture  iv  =  r^/r^  =  0.39).  Also  shown  in  Fig.  III-2,  is  an  example  of  the  instru¬ 
mental  MTF  degraded  by  the  atmosphere.  The  case  shown  is  for  a  5.4-km  horizontal  path  and 
was  taken  at  13:10  EDT.  The  results  of  a  more  comprehensive  set  of  measurements  on  the 
atmospheric  MTF,  made  in  conjunction  with  10.6-|jLm  propagation  measurements  also  described 
herein,  will  be  presented  in  the  next  Optics  Research  Report. 

J.  Z.  Holtz 
T.  J.  Gilmartin 
L.  J,  Prugnarola 


2.  MTF  Electronics 

The  electronics  used  with  the  MTF  instrument  consists  of  a  square-wave  generator  and  a 
signal  processor. 

The  square-wave  generator  consists  of  a  light  source  and  photodetector  in  an  integral  hous¬ 
ing  with  a  beam  splitter  and  lens  such  that  the  light  is  focused  on  the  rotating  reticle  and  the  re¬ 
turn  beam  is  focused  on  the  photodetector.  Following  the  photofet  detector  are  an  amplifier  and 
an  emitter  follower  which  provides  a  low  output  impedance.  The  square-wave  generator  pro¬ 
vides  the  reference  signal  when  a  lock-in  amplifier  is  used  to  process  the  data. 

The  signal  processor  provides  an  output  proportional  to  the  percentage  modulation  im¬ 
pressed  on  the  incoming  beam  by  the  rotating  reticle.  It  includes:  * 

(a)  Photomultiplier 

(b)  Photomultiplier  buffer  4 

(c)  Low-pass  active  filter  before  the  signal  strength  meter 

(d)  AGC  loop  amplifier  and  filter,  which  control  the  voltage  applied  to 
the  photomultiplier  and  thus  its  gain 
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(e)  Active  bandpass  filter  tuned  to  the  frequency  produced  by  the  bar 
pattern  of  the  rotating  reticle 

(f)  Full- wave  rectifier  and  filter,  which  convert  the  output  of  the  bandpass 
filter  into  a  DC  voltage  proportional  to  the  percentage  modulation. 


ROTATING 

RETICLE 


WAVE 


Fig.  111-3.  MTF  control  and  signal  processing  electronics. 

Outputs  are  available  at  various  points  for  recording- or  for  connecting  other  measuring  instru¬ 
ments.  Figure  III- 3  is  a  block  diagram  of  the  electronics.  Zimmerman 


B.  AUTOMATIC  RANGE  TRACKER 

Detailed  plans  for  an  automatic  range  tracker  for  the  Firepond  laser  radar  installation  were 
devised.  Implementation  of  the  plans  has  been  cancelled  because  of  planned  program  and  system 
changes  at  Firepond. 

An  automatic  switch  designed  to  interchange  the  tuning  relationships  of  the  local  oscillator 
(LO)  to  the  master  oscillator  (MO)  as  controlled  by  the  delta  f  loop  servo  has  been  constructed 
and  is  ready  to  be  installed  and  tested.  It  will  switch  the  output  frequency  of  the  LO  from  5  MHz 
below  the  MO  frequency  to  5  MHz  above  the  MO  frequency  when  the  Doppler  frequency  goes  to 
zero.  In  this  manner,  the  signal  frequency  at  the  detector  mixer  output  will  be  prevented  from 
going  to  zero  regardless  of  whether  the  moving  target  is  on  the  incoming  or  outgoing  leg  of  its 
pass.  Figure  III-4  is  a  block  diagram  of  the  switching  scheme. 

Design  and  construction  of  the  remote  control  circuits  to  be  used  with  the  AIL  monopulse 
tracking  receiver  have  been  completed.  The  circuits  are  ready  to  be  installed  and  tested  at 
Firepond  when  the  receiver  arrives.  Zimmerman 
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Fig. 111-4.  Automatic  Doppler  switchover. 

C.  ATMOSPHERIC  OPTICAL  COHERENCE  MEASUREMENTS  AT  10.6  AND  0.63  jLim 

The  performance  of  large  aperture  optics  (D  >  Im)  at  10.6  pm  in  projection  and  heterodyne 
reception  has  previously  been  estimated  by  extrapolating  from  measurements  made  at  shorter 
wavelengths.  Since  large  and  complex  systems  are  being  contemplated  at  10.6  pm,  it  is  of  great 
practical  interest  to  measure  directly  the  effect  which  atmospheric  refractive  inhomogeneities 
will  have  on  such  systems.  To  this  end,  a  48~inch-diameter  (1.2  m)  Cassegrainian  telescope 
was  focused  at  a  range  of  5.4 km^  where  the  10.6-pm  spot  fell  on  a  diffuse  metal  surface,  was 
scanned  with  an  IR  imaging  camera  and  recorded  for  later  computer  analysis.  Simultaneously, 
the  optical  quality  of  parallel  paths  1  and  2  m  distant  from  the  primary  path  were  monitored  at 
0.63  pm  with  a  MTF  instrument  and  with  a  two-point,  white  light  source  resolution  instrument. 

The  optical  paths  were  horizontal  with  an  average  altitude  of  68  m. 

Data  have  now  been  taken  under  a  variety  of  meteorological  conditions  and  over  several 
diurnal  cycles.  However,  while  most  of  the  data  from  the  visible  instruments  have  been  analyzed, 
the  computer  analysis  of  the  10.6-pm  data  will  not  be  completed  for  several  weeks.  For  this 
reason,  a  detailed  description  of  the  experiment  or  of  the  results  will  not  be  presented  at  this 
time,  but  some  qualitative  observations  will  be  made. 

Defining  the  atmospheric  coherence  diameter  r^  to  be  the  diameter  which,  when  used  in  a. 
vacuum,  allows  the  same  resolution  as  does  the  atmosphere  when  not  aperture  limited,  the  at¬ 
mospheric  r  was  measured  by  the  MTF  instrument  to  range  from  about  1  to  20  cm  at  0.63  pm. 
Since  r  is  predicted  to  scale  as  \  ^  ,  r^(10.6  pm)  is  expected  to  equal  29.4  r^(0.63  pm)  and  the 
r  (10.6  pm),  implied  by  the  visible  measurements,  ranges  from  30  cm  to  6  m.  Indeed,  the  image 
projected  by  the  1.2-m  aperture  at  10.6  pm  was  observed  to  be  almost  perfectly  diffraction  lim¬ 
ited  occasionally,  near  diffraction  limited  most  often,  and  badly  degraded  occasionally.  A 
measurement-by- measurement  correlation  of  coherence  data  taken  at  the  two  wavelengths  is 
presently  being  made. 

The  10.6-pm  spot  wander  was  observed  to  range  from  a  fraction  of  the  diffraction-limited 
spot  radius  to  several  spot  radii  at  extremes  under  typical  wind  and  turbulence  conditions.  The 
correlation  time  for  spot  wander  and  image  distortion  was  typically  a  fraction  of  a  second. 
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A  more  detailed  and  comprehensive  description  of  the  results  of  these  propagation  measure¬ 


ments  will  be  given  in  the  next  Optics  Research  Report. 


T,  J.  Gilmartin 
J.  Z.  Holtz 
R.  N.  Capes 


D.  InSb  ISOLATOR 

We  have  previously  reported  that  when  10. 6 -pm  radiation  is  incident  on  an  InSb  sample  that 
is  rigidly  attached  to  a  heat  sink,  beam  distortion  occurs.  The  stress  caused  by  local  heating 
will  not  be  uniform  in  this  case  and  the  resultant  unsymmetrical  change  in  the  index  of  refraction 
will  cause  the  distortion. 

A  new  holder  has  been  designed  which  allows  the  sample  to  expand  uniformly.  The  holder 
is  shown  in  Fig.  III-5.  Tests  using  a  12-W  beam  incident  on  a  4-mm  piece  of  InSb  held  by  the 
new  holder  show  none  of  the  distortion  previously  observed. 


Fig.lll-5.  InSb  holder. 


A  profilometer  was  constructed  which  scans  the  beam  at  if-inch/min.  (25  thousandths/sec) 
through  a  20  to  30  thousandths  aperture.  The  time  constant  of  the  detector  is  10  sec.  Results 
from  the  profilometer  show  no  gross  differences  between  the  beam  directly  from  the  laser  and 
one  that  has  passed  through  an  InSb  sample. 

An  isolator  was  constructed  using  Faraday  free  carrier  rotation.  Pictures  of  the  system 
are  shown  in  Figs.  III-6  and  -7.  The  Dewar  and  a  polarizer  are  shown  in  Figs.  III-8  and  -9, 
res'pectively. 

The  isolator  has  been  tested  in  the  laboratory.  Characteristics  of  the  10.6-pm  isolator  are: 


Material 

Concentration 

Thickness 

Verdet  constant 
Field  for  45  °  rotation 
Operating  temperature 


n  —  InSb 

1.0  X  10^^  cm"^ 

4  mm 

1.1  X  10"^deg/cm-G 

10.2kG 

77°K 
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Fig.lll“8.  Exploded  view  of  Dewar:  on  left, 
the  outer  Dewar  with  Ge  A-R  coated  windows; 
on  right,  the  inner  liquid  nitrogen  reservoir  with 
InSb  attached.  Liquid  nitrogen  flows  to  the  cop¬ 
per  plug  visible  [ust  below  the  InSb  holder. 
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Fig.  1 11-9.  Exploded  view  of  Ge  brewster  and  polarizer. 
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Insertion  loss 
Isolation 

Total  reverse  loss 

Maximum  forward 
throughput  power 

Safety  factor 


2  ±  O.BdB 
20  ±  O.SdB 
22  ±  idB 
SOW  cw 

2.9  (calculated) 


The  safety  factor  is  the  ratio  of  the  radius  of  any  aperture  to  the  beam  waist  at  the  worst  point 
in  the  optical  path.  In  this  case  it  occurs  at  the  point  where  the  pole  faces  and  the  magnet  coils 
meet. 

Isolation  seems  to  be  limited  by  non -uniformity  in  the  magnetic  field  across  the  ^-inch 
aperture  in  the  pole  piece.  Maximum  isolation  based  on  the  limitation  of  the  Ge  brewster  angle 
polarizers  is  26 dB.  [There  is  a  23  percent  transmission  of  perpendicular  polarization  through 
each  surface.  There  are  4  surfaces,  thus  the  isolation  is  0.23"^.]  Notice  we  are  6dB  below  this 
level  indicating  that  this  is  not  the  limiting  factor. 

Some  improvement  might  be  possible  by  using  a  magnet  with  larger  pole  faces.  Larger 
pole  faces  generally  reduce  variation  of  field  across  a  :^-inch-diameter  hole  in  the  pole  faces. 

An  investigation  of  the  power  handling  capability  of  the  isolator  when  300-nsec  pulses  are 
incident  will  be  undertaken.  The  effect  of  InSb  on  the  waveform  of  the  pulse  will  be  examined. 


L.  R.  Tomas etta 
W.  E.  Bicknell 
D.  H.  Bates 


i 


i 


92 


UNCLASSIFIED 
Security  Classification 


DOCUMENT  CONTROL  DATA  -  R&D 

(Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  is  classified) 


1.  ORIGINATING  ACTIVITY  (Corporate  author) 

Lincoln  Laboratory,  M.I.T. 

3.  REPORT  TITLE 

Optics  Research 


4.  DESCRIPTIVE  NOTES  (Type  of  report  and  inclusive  dates) 

Semiannual  Report  —  1  January  through  30  June  1972 

5.  AUTHOR(S)  CLast  name,  first  name,  initial) 

Rediker,  Robert  H. 


6.  REPORT  DATE 

30  June  1972 


8a.  CONTRACT  OR  GRANT  NO.  F 19628-70-C-0230 
fc.  PROJECT  NO.  ARPA  Order  600 


10.  AVAILABILITY/LIMITATION  NOTICES 


Approved  for  public  release;  distribution  unlimited. 


2a.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

2b.  GROUP 

None 


7a.  TOTAL  NO.  OF  PAGES 

110 

76.  NO.  OF  REFS 

69 

9a.  ORIGINATOR'S  REPORT  I 

slUMBER(S) 

Optics  Research  (1972:1) 

96.  OTHER  REPORT  NO(S)  (Any  other  numbers  that  may  be 

assigned  this  report) 

ESD-TR-72-195 

11.  SUPPLEMENTARY  NOTES 

None 


12.  SPONSORING  MILITARY  ACTIVITY 


Advanced  Research  Projects  Agency, 
Department  of  Defense 


13.  ABSTRACT 


This  report  covers  work  of  the  Optics  Division  at  Lincoln  Laboratory  for 
the  period  1  January  through  30  June  1972.  The  topics  covered  are  laser  tech¬ 
nology  and  propagation,  optical  measurements  and  instrumentation,  and  laser 
radar.  Additional  information  on  the  optics  program  may  be  found  in  the  semi¬ 
annual  technical  summary  reports  to  the  Advanced  Research  Projects  Agency. 


14.  KEY  WORDS 


optics 

laser  technology 
thermal  blooming 


optical  devices 
optical  systems 
KC-135 


laser  radar  imaging  system 
air  pollution  control 


UNCLASSIFIED 
Security  Classification 


